First and Second Order Equations for Piezoelectric 
Crystals Expressed in Tensor Form 

By W. P. MASON 

Introduction 

AEOLOTROPIC substances have been used for a wide variety of elastic 
piezoelectric, dielectric, pyroelectric, temperature expansive, piezo- 
optic and electro-optic effects. While most of these effects may be found 
treated in various publications there does not appear to be any integrated 
treatment of them by the tensor method which greatly simplifies the method 
of writing and manipulating the relations between fundamental quantities. 
Other short hand methods such as the matrix method can also be used for 
all the linear effects, but for second order effects involving tensors higher 
than rank four, tensor methods are essential. Accordingly, it is the purpose 
of this paper to present such a derivation. The notation used is that agreed 
upon by a committee of piezoelectric experts under the auspices of the Insti- 
tute of Radio Engineers. 

In the first part the definition of stress and strain are given and their inter- 
relation, the generalized Hookes law is discussed. The modifications caused 
by adiabatic conditions are considered. When electric fields, stresses, and 
temperature changes are applied, there are nine first order effects each of 
which requires a tensor to express the resulting constants. The effects are 
the elastic effect, the direct and inverse piezoelectric effects, the temperature 
expansion effect, the dielectric effect, the pyroelectric effect, the heat of 
deformation, the electrocaloric effect, and the specific heat. There are 
three relations between these nine effects. Making use of the tensor trans- 
formation of axes, the results of the symmetries existing for the 32 types of 
crystals are investigated and the possible constants are derived for these 
nine effects. 

Methods are discussed for measuring these properties for all 32 crystal 
classes. By measuring the constants of a specified number of oriented cuts 
for each crystal class, vibrating in longitudinal and shear modes, all of the 
elastic, dielectric and piezoelectric constants can be obtained. Methods 
for calculating the properties of the oriented cuts are given and for deriving 
the fundamental constants from these measurements. 

1 For example Voigt, "Lehrbuch der Kristall Physik," B. Teubner, 1910; Wooster, 
"Crystal Physics," Cambridge Press, 1938; Cady "Piezo-electricity" McGraw Hill, 1946. 
3 The matrix method is well described by W. L. Bond "The Mathematics of the Physical 
Properties of Crystals," B. S. T. J., Vol. 22, pp. 1-72, 1943. 
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Second order effects are also considered. These effects (neglecting second 
order temperature effects) are elastic constants whose values depend on 
the applied stress and the electric displacement, the electrostrictive effect, 
piezoelectric constants that depend on the applied stress, the piezo-optical 
effect and the electro-optical effect. These second order equations can 
also be used to discuss the changes that occur in ferroelectric type crystals 
such as Rochelle Salt, for which between the temperature of — 18°C. and 
+ 24°C, a spontaneous polarization occurs along one direction in the crystal. 
This spontaneous polarization gives rise to a first order piezoelectric deforma- 
tion and to second order electrostrictive effects. It produces changes in 
the elastic constants, the piezoelectric constants and the dielectric constants. 
Some measurements have been made for Rochelle Salt evaluating these 
second order constants. 

Mueller in his theory of Rochelle Salt considers that the crystal changes 
from an orthorhombic crystal to a monoclinic crystal when it becomes 
spontaneously polarized. An alternate view developed here is that all of 
the new constants created by the spontaneous polarization are the result of 
second order effects in the orthorhombic crystal. As shown in section 7 
these produce new constants proportional to the square of the spontaneous 
polarization which are the ones existing in a monoclinic crystal. On this 
view "morphic" effects are second order effects produced by the spontaneous 
polarization. 

1. Stress and Strain Relations in Aeolotropic Crystals 

I.I. Specification of Stress 

The stresses exerted on any elementary cube of material with its edges 
along the three rectangular axes X, Y and Z can be specified by considering 
the stresses on each face of the cube illustrated by Fig. 1. The total stress 
acting on the face ABCD normal to the X axis can be represented by a 
resultant force R, with its center of application at the center of the face, 
plus a couple which takes account of the variation of the stress across the 
face. The force R is directed outward, since a stress is considered posi- 
tive if it exerts a tension. As the face is shrunk in size, the force R will be 
proportional to the area of the face, while the couple will vary as the cube of 
the dimension. Hence in the limit the couple can be neglected with respect 
to the force R. The stress (force per unit area) due to R can be resolved 
into three components along the three axes to which we give the designation 

■* "2 > "* V*2 i * ZX 2 • (1 ) 

Here the iirst letter designates the direction of the stress component and the 
second letter .r 2 denotes the second face of the cube normal to the X axis. 
Similarly for the first X face OEFG, the stress resultant can be resolved 
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into the components T XXl , T yXl , T ZTJ , which are oppositely directed to 
those of the second face. The remaining stress components on the other 
four faces have the designation 



(2) 




Fig. 1. — Cube showing method for specifying stresses. 

The resultant force in the X direction is obtained by summing all the forces 
with components in the X direction or 

F x = (T XXl - T XXi ) dydz + (T m - T m ) dxdz + (T xgJ - T xtl ) dxdy. (3) 



But 



dT x 



T t x, — ~T X Xi + —t— dx; lxy 2 — 1*V\ 



+ if dy] . r -" ~ T m +""'!* 



djj, 

dz 



(4) 



and equation (3) can be written in the form 



F. = 



--( 



dT xx , dTx V , dT x 



dx dy dz. 



dx dy dz 

Similarly the resultant forces in the other directions are 

(dTyX , dTyy . dJ\ 

\ dx dy dz 



(5) 



p = -( xj*? + zis + ^ j dx dy dz 



■) 



*--(t?.+i? +*)*** 



(6) 
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We call the components 



T, 



T X2 




Tu, 


Tn, 


T n 


Ty, 


— 


Tn, 


Tvii 


T» 


T„ 




Tn, 


T&, 


T 33 



(7) 



the stress components exerted on the elementary cube which tend to deform 
it. The rate of change of these stresses determines the resultant force on 
the cube. The second form of (7) is commonly used when the stresses are 
considered as a second rank tensor. 




Fig. 2.— Shearing stresses exerted on a cube. 

It can be shown that there is a relation between 3 pairs of these compo- 
nents, namely 



' xv — 1 v 



r«. = r, 



I yt — I m • 



(8) 



To show this consider Fig. 2 which shows the stresses tending to rotate the 
elementary cube about the Z axis. The stresses T yxj and T yXl tend to rotate 
the cube about the Z axis by producing the couple 



T yX dx dy dz 



(9) 



The stresses T ZVl and T XVl produce a couple tending to cause a rotation in 
the opposite direction so that 



5 {T yx — Txy) dx dy dz = couple = Iib z 



(10) 



is the total couple tending to produce a rotation around the Z axis. 
But from dynamics, it is known that this couple is equal to the product of 
the moment of inertia of the section times the angular acceleration. This 
moment of inertia of the section is proportional to the fourth power of the 
cube edge and the angular acceleration is finite. Hence as the cube edge 
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* *V 1 


T xt 


\Tu, 


Tn, 


T u 




Ti, 


n, 


r. 


■* VV > 


Ty X 


= \Tu, 


T22 > 


r 23 
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T; 


Tu 


T A 


■* 1/2 J 


T„ 


| ^13 >' 


2^23 > 


T33 




T tl 


Ta, 


Ti 



approaches zero, the right hand side of (10) is one order smaller than the 
left hand side and hence 

Ty^Tty. (11) 

The same argument applies to the other terms. Hence the stress com- 
ponents of (7) can be written in the symmetrical form 



(12) 



The last form is a short hand method for reducing the number of indices 
in the stress tensor. The reduced indices 1 to 6, correspond to the tensor 
indices if we replace 

llbyl; 22 by 2; 33 by 3; 23 by 4; 13 by 5; 12 by 6. 

This last method is the most common way for writing the stresses. 

1.2 Strain Components 

The types of strain present in a body can be specified by considering two 
points P and Q of a medium, and calculating their separation in the strained 
condition. Let us consider the point P at the origin of coordinates and the 
point Q having the coordinates x, y and z as shown by Fig. 3. Upon strain- 




Fig. 3. — Change in length and position of a line due to strain in a solid body. 

ing the body, the points change to the positions P', Q'. In order to specify 
the strains, we have to calculate the difference in length after straining, or 
have to evaluate the distance P'Q'-P Q. After the material has stretched 
the point P' will have the coordinates & , in , ft , while Q' will have the 
coordinates x + £2 ; y + V2 ; z + f 2 . But the displacement is a continuous 
function of the coordinates x, y and z so that we have 
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Similarly 

, drj . dr] . dr] 
™ = '' + F x X+ 3y y + lTz Z 

Hence subtracting the two lengths, we find that the increases in separation 
in the three directions are 

*-•? + *£ + ■? (14) 

y dx dy dz 

ar , ar , ar 

The net elongation of the line in the .r direction is a — and the elongation 

per unit length is — which is defined as the linear strain in the x direction. 

c dx 

We have therefore that the linear strains in the x, y and 3 directions are 

S «« ; &-*>; *-£. (15) 

dx dy dz 

The remaining strain coefficients are usually defined as 

dy dz dz dx dx dy 

and the rotation coefficients by the equations 

u * = % ~ 5 ; w " = a~ 3 " te ' "* ~ & ay" (l7) 

Hence the relative displacement of any two points can be expressed as 

K = x (^|^) + > & + * (^T^) (18) 

,. = . (&Z*) + y (* + ») + S5l 
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which represents the most general type of displacement that the line P Q 
can undergo. 

As discussed in section 4 the definition of the shearing strains given by- 
equation (16) does not allow them to be represented as part of a tensor. 
If however we defined the shearing strains as 



2S 23 = S< = ^ + £) ; 25,3 - S b 



they can be expressed in the form of a symmetrical tensor 



(19) 
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(20) 



For an element suffering a shearing strain S 6 = 2Sn only, the displace- 
ment along x is proportional to y, while the displacement along y is propor- 
tional to the x dimension. A cubic element of volume will be strained into 
a rhombic form, as shown by Fig. 4, and the cosine of the resulting angle 6 




Fig. 4. — Distortion due to a shearing strain. 

measures the shearing deformation. For an element suffering a rotation 
co 2 only, the displacement along x is proportional to y and in the negative 
y direction, while the displacement along y is in the positive x direction. 
Hence a rectangle has the displacement shown by Fig. 5, which is a pure 
rotation of the body without change of form, about the z axis. For any 
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body in equilibrium or in nonrotational vibration, the co's can be set equal 
to zero. 

The total potential energy stored in a general distortion can be calculated 
as the sum of the energies due to the distortion of the various modes. For 

example in expanding the cube in the x direction by an amount — dx = 

Si dx, the work done is the force times the displacement. The force wil 




Fig. 5. — A rotationbf a solid body. 

be the force 7\ and will be T\ dy dz. Hence the potential energy stored in 
this distortion is 



T\ dSi dx dy dz 
For a shearing stress T 6 of the type shown by Fig. 4 the displacement 



dStdx 



J C 1 J 

times the force T 6 dy dz and the displacement 6 times the force T t dx dz 
equals the stored energy or 

APEz = | (dS 6 T 6 + dSaT 6 ) dx dy dz = dSsT 6 dx dy dz. 
Hence for all modes of motion the stored potential ener gy is equal to 
APE = [Ti dSi + T 2 dS 2 + T 3 dS 3 + T 4 dS 4 + T & dS B 



+ Tt dSi] dx dy dz. 



(21) 



1.3 Generalized Hooke's Law 



Having specified stresses and strains, we next consider the relationship 
between them. For small displacements, it is a consequence of Hooke's 
Law that the stresses are proportional to the strains. For the most un- 
symmetrical medium, this proportionality can be written in the form 



(22) 
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T\ = C\\S\ + CjsSi -\- C13S3 -\- C1404 "T C15VJ5 + C16O6 

7*2 =• C21S1 -\~ C22S2 ~\~ C23S3 -\- C04S4 + £26^6 I C2frS*6 

7^3 = C31.S1 + C32S2 ~f" C33S3 + C3iS4 + £36^6 "T" ^B^G 

7*4 = C41^1 + C42iS2 -f- C4353 + £44*!) 4 + C45O5 -f" C46»->6 

^b = C51S1 ~T" CssSi + CB3O3 + C54O4 T C55O5 + C56O6 

7*6 = CgK>1 H~ C62S2 4" C63O3 + C64O4 T" C&&S& T" ^66^6 

where Cu for example is an elastic constant expressing the proportionality 
between the Si strain and the T\ stress in the absence of any other strains. 
It can be shown that the law of conservation of energy, it is a necessary 
consequence that 

C12 = C21 and in general dj = ca . (23) 

This reduces the number of independent elastic constants for the most 
unsymmetrical medium to 21. As shown in a later section, any symmetry 
existing in the crystal will reduce the possible number of elastic constants 
and simplify the stress strain relationship of equation (22). 

Introducing the values of the stresses from (22) in the expression for the 
potential energy (21), this can be written in the form 

2PE = CuSl + 2c 12 5i5 2 + 2c 13 S!S s + 2c li S 1 S i + 2c 16 S 1 S 6 + 2a 6 S 1 S 6 

-+" C22S2 "T" 2C23S2S3 + 2C24S2S4 + 2C2bS2iSb "T" 2C26*S , 206 

~r C3303 -j- 2C34O304 -j- 2^350365 -f- 2C36O3O6 

+ C44.S4 -+- 2c455V>5 + 2^460406 (24) 

+ C6frS5 + 2c5o5556 
+ ^66^6 • 

The relations (22) thus can be obtained by differentiating the potential 
energy according to the relation 

T - dPE - • T - dPE (?K\ 

doi doe 

It is sometimes advantageous to express the strains in terms of the stresses. 
This can be done by solving the equations (22) simultaneously for the 
strains resulting in the equations 
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Si = SnTi + 5i 2 r 2 + S13T3 + SuTi + Sib^b + 5i 6 r 6 

52 = SnTi + 522^2 + S23T3 + SuTi + S2 & T h + $26^6 

53 = S31T1 + 5 3 2^2 + ^33^3 + $34^4 + fts^B + S^Tt 

54 = SnTi + 542^2 + S43T3 + 5 44^4 + $45^5 + SitTf 
Sf, = S b \Ti + 562^2 + S^Ta + ^54^4 + S 66 T 6 + S 66 T 6 
Sq = S6lT\ + 562^2 + 563^3 + 564^4 + 566^6 + ^66^6 



(26) 



where 



i+j 



Sii = 



(-d^a:-,- 



(27) 



for which A c is the determinant of the c, ; - terms of (28) and A<j the minor 
obtained by suppressing the *th and^'th column 



A c = 



C\\ C\3 C13 • Cu C15 fl6 

fl2 C22 ^23 ^24 C-2b C25 

Cl3 Cn ^33 C34 C35 C36 

Cu C24 C34 C44 C45 f 46 

Cl6 C2a ^35 C45 C55 C56 

Cl6 ^26 ^36 C46 ^66 ^66 



(28) 



Since Ctj = Cji it follows that s,-,- = 5,-,-. The potential energy can be 
expressed in the form. 

2PE = s n Tl + TsnTiTi + IhFiTi + 2suTiTi + 2* 1B r 1 r 5 + Is^T* 

+ ^22^2 + 2S23T2T3 + 2s u TzTi + 2si{TzTi + 2si f ,T?,T( l 

+ *„7l + 2534r 3 r 4 + 2s 35 r 3 rB + 2*8 fl r s r« 

+ 544^4 + 2545r 4 r 5 + 2SwTiT 6 

+ SbaTl + 2556^5^6 

+ $66-' 6 • 

The relations (26) can then be derived from expressions of the type 



(29) 



Si = 



dPE 



S& = 



dPE 



(30) 



dTi ' ' dT 6 

1.4 Isothermal and Adiabatic Elastic Constants 

We have so far considered only the elastic relations that can be measured 
statically at a constant temperature. The elastic constants are then the 
isothermal constants. For a rapidly vibrating body, however, there is no 
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chance for heat to equalize and consequently the elastic constants operative 
are the adiabatic constants determined by the fact that no heat is added 
or subtracted from any elemental volume. For gases there is a marked 
difference between the adiabatic and the isothermal constants, but for 
piezoelectric crystals the difference is small and can usually be neglected. 

To investigate the relation existing we can write from the first and second 
laws of thermodynamics, the relations 

dU = [7i dSi + T 2 dS 2 + T 3 dS 3 

(31) 

+ r 4 dSt + t 6 ds 6 + r 6 dSs] + e<h 

which expresses the fact that the change in the total energy U is equal to 
the change in the potential energy plus the added heat energy dQ = Q da 
where is the temperature and a the entropy. Developing the strains and 
entropy in terms of the partial differentials of the stresses and temperature, 
we have 

dS, = ffr 1 «■■ + fl 1 4T, + §'dT: 

Oil Oil Ol 3 

+ W, dTt + if. dTi + if. dT ' + ae de 



dSfi Trr , . dSs 1rrl i dSt 



dS< = £*dTi + ~'tfTa + ^dT, 

Oil olz oli 



+ Wr i dTl + m dTi + Bf t dT ' + m dB 



(32) 



The partial derivatives of the strains with regard to the stresses are readily 
seen to be the isothermal elastic compliances. The partial derivatives of 
the strains by the temperatures are the six temperature coefficients of ex- 
pansion, or 

dSi . dS 6 . . 

ae""*' ""' de = " 6 - (33) 

To evaluate the partial derivatives of the entropy with respect to the 
stresses we make use of the fact that U is a perfect differential so that 

dSi da dSt da f%*\ 

de = dT 1 = ai; ""' de = dT t = a6 - (34) 
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Finally multiplying through the last of equation (32) by 9 we can write 

them as 

5, = tfiJTi + *? 2 r 2 + 5? 8 r 3 + slT, + sf 6 T b + s? e T 6 + ax <JB 

(35) 
St = sftTx + s%T 2 + s? 6 r 3 + J?er 4 + fftTi + Aft + a 6 <*9 

JQ = 6 Jo- = e[air a + a 2 r 2 + a 3 r 3 + cuT* + a B r 6 + a 6 T 6 ] + P C„ </9 

since 9 — is the total heat capacity of the unit volume at constant stress, 

6*9 
which is equal to pC p , where p is the density and C p the heat capacity at 
constant stress per gram of the material. 

To get the adiabatic elastic constants which correspond to no heat loss 
from the element, or dQ = 0, dO can be eliminated from (35) giving 

Si = siiTi + s'nTi + sl 3 T 3 + s' u T 4 + s{,,Tt + s' u T t + WpC p ) dQ 
(36) 

S 6 = s'uTi + s' it T a + s"z 6 T 3 + s' M T 4 + slsTs + s' M T t + (cf/pCp) dQ 

where 

, e «j a, 9 ,.,->. 

pL p 
For example for quartz, the expansion coefficients are 

on = 14.3 X 10- e / o C; a, = 14.3 X KT 6 /°C; «s - 7.8 X l(r7°C; 

i 

Oii = «5 = a 6 = 

The density and specific heat at constant pressure are 

p = 2.65 grams/cm 8 ; C p = 7.37 X 10 6 ergs/cm 3 . 
Hence the only constants that differ for adiabatic and isothermal values are 

Su = s 22 ; S12; S13; S33 ■ 
Taking these values as 

s'n = 127.9 X 10" 14 cm 2 /dyne; s' a = -15.35 X 10 -14 ; 
s'u = 11.0 X 1(T 14 ; s' K = 95.6 X 10 -14 . 
We find that the corresponding isothermal values are 

sti = 128.2 X 1(T 14 ; s? 2 = -15.04 X 10 -14 ; 
5? 3 = 10.83 X 10" 14 ; s? 8 = 95.7 X 10 -14 cm 2 /dyne 

3 See "Quartz Crystal Applications" Bell System Technical Journal, Vol. XXII. 
No. 2, July 1943, W. P. Mason. 
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at 25°C. or 298° absolute. These differences are probably smaller than 
the accuracy of the measured constants. 

If we express the stresses in terms of the strains by solving equation (35) 
simultaneously, we find for the stresses 

Ti = cnSi + CnSi + c u S 3 + c u Si + ci 6 S 6 + Ci 6 Si — Xi dQ 
(38) 

T 6 = Cieoi -p C26-J2 I C36O3 -\~ C4804 T Cb&S b I CgeOa — Afl dQ 

where 

\ -6 1 61 61 61 0| 9 

Ai = aidi -J- 02^12 t a 3 Ci3 -+- a.iCu + a b Ci b + a 6 Ci 6 



\ e 1 e 1 o 1 61 Oi e 

X ti = «lCi6 + 0^26 + a 3 C 36 + 04C 4 6 + OC b C b6 + (X 6 C 66 . 

The X's represent the temperature coefficients of stress when all the strains 
are zero. The negative sign indicates that a negative stress (a compression) 
has to be applied to keep the strains zero. If we substitute equations (38) 
in the last of equations (35), the relation between increments of heat and 
temperature, we have 

dQ = Qda = 0[aiSi + x 2 s 2 + a 3 s 3 + X4S4 + \ b S 6 + X 6 S 6 ] 

(39) 
+ [pCp — 0(aiXi + 0:2X2 + a 3 X 3 + 04X4 + "5X5 + a 6 \^)]dQ. 

If we set the strains equal to zero, the size of the element does not change, 
and hence the ratio between dQ and dQ should equal p times the specific 
heat at constant volume C v . We have therefore the relation 

p[C p — C v ] = 6[q!iXi + 0:2X2 + 0:3X3 + 0:4X4 + 0:5X5 + qgXe]' (W) 

The relation between the adiabatic and isothermal elastic constants c,y 
thus becomes 

Cu = c ti + —pi— . (41) 

pL v 

Since the difference between the adiabatic and isothermal constants is so 
small, no differentiation will be made between them in the following sections. 

2. Expression for The Elastic, Piezoelectric, Pyroelectric and 
Dielectric Relations of a Piezoelectric Crystal 

When a crystal is piezoelectric, a potential energy is stored in the crystal 
when a voltage is applied to the crystal. Hence the energy expressions of 
(31) requires additional terms to represent the increment of energy dU. 
If we employ CGS units which have so far been most widely used, as applied 
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lo piezoelectric crystals, the energy stored in any unit volume of the crystal is 
dU = Ti dSi + T 2 dS 2 + T 3 dS 3 + T 4 dS* + T b dS b + T 6 dS 6 

. _ dD, , __ dD 2 , -, </D 3 , Q , < 42) 

47T 47T 47T 

where £1 , £2 and £3 are the components of the field existing in the crystal 
and A , D 2 and D 3 the components of the electric displacement. In order 
to avoid using the factor l/4ir we make the substitution 

£ - 5. (43) 

47T 

The normal component of 8 at any bounding surface is S the surface charge. 
On the other hand if we employ the MKS systems of units the energy of 
any component is given by E n dD n directly and in the following formulation 5 
can be replaced by D. 

There are two logical methods of writing the elastic, piezoelectric, pyro- 
electric and dielectric relations. One considers the independent variables 
as the stresses, fields, and temperature, and the dependent variables as the 
strains, displacements and entropy. The other system considers the strains, 
displacements and entropy as the fundamental independent variables and 
the stresses, fields, and temperature as the independent variables. The 
first system appears to be more fundamental for ferroelectric types of 
crystals. 

If we develop the stresses, fields, and temperature in terms of their partial 
derivatives, we can write 

*-£) « + sP) <* + £) <* + §).* 

dS\/D,a 002/ D.a Ob 3 /D.<r OO4/ D.a 

)) *> + %) « + £) * + £) « 

b/D,<r Ob 6 /D,a OOl/S.a O0 2 /8,o 



T <~ as 



003/ S,a Off /S.D 



(44 A) 



l/D,o Ob*/ D.a abz/D,a OO4/ D.ff 

) *. + £•) «+£•) *+S) « 

/D.o Obf./U.a Odi/s.a Od 2 /S,a 

) d!k+ a -p) da 

/S,a oa /S,D 



+ t> '«-- 



'~ d8 3 
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£, = £,= °£) iS, + H) iS, + §f) iS, + *f ) dS. 

OOl/D.a O02/D,* OOZ/D.0 OOt/D.a 

4- —±\ dS 4- — • I dS 4- ^M <f5 4- — l i rf6 

dSf,/D,a dSs/D.ir d8i/B,<r d&2 / B,o 



+ «£\ i S , + 3 ^) i. 

T /s,n 



55a /a,<r d 



OOl/D.a 002/0.ff OOa/D.a OOi/D.a 

') dS* + d £) iS. + ^f) <% + 

5/D.o OOQ/D.o OOl /B,o 

*>-£) "'+£) *■+£) «+k) « 

OOl/D,t a02/D,c OOa/D.a OOi/D.o 



Odf,/D,o OOQ/D.o OOl /B,o 002/3,<r 

+ «) * + «*) & < 44B > 

003/ B.c Off /B.D 



+ ") «+«) «,+£) *+£) * 

UO5/ D,o Od6/D,o O0i/8,o O02 /S.b 

60 s 



*) * + *) «,. 

003/ B,a Off/S,D 



38, 

The subscripts under the partial derivatives indicate the quantities kept 
constant. A subscript Z> indicates that the electric induction is held 
constant, a subscript a indicates that the entropy is held constant, while a 
subscript S indicates that the strains are held constant. 

Examining the first equation, we see that the partial derivatives of the 
stress T\ by the strains are the elastic constants c,-,- which determine the 
ratios between the stress Ti and the appropriate strain with all other strains 
equal to zero. To indicate the conditions for the partial derivatives, the 
superscripts D and a are given to the elastic constants and they are written 
Cif. The partial derivatives of the stresses by 8 = D/4r are the piezo- 
electric constants /?,-, which measure the increases in stress necessary to 
hold the crystal free from strain in the presence of a displacement. Since 
if the crystal tends to expand on the application of a displacement, the 
stress to keep it from expanding has to be a compression or negative stress, 
the negative sign is given to the k'j constants. As the only meaning of 
the h constants is obtained by measuring the ratio of the stress to 8 = D/lir 
at constant strains, no superscript S is added. However there is a difference 
between isothermal and adiabatic piezoelectric constants in general, so 
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that these piezoelectric constants are written /;*•;. Finally the last partial 
derivatives of the stresses by the entropy a can be written 

d -^) d ° = kir) Qd ° = iir) d Q = ~y^ d Q («) 

da /s.d 6 da /s,d W oar /b,d 

where dQ is the added heat. We designate 1/6 times the partial derivative 
as — 7n° and note that it determines the negative stress (compression) 
necessary to put on the crystal to keep it from expanding when an increment 
of heat dQ is added to the crystal. The electric displacement is held 
constant and hence the superscripts S, and D are used. The first six equa- 
tions then can be written in the form 

T n = CnlSi + C„£S2 + Cn'zSz + C^Si + C„'i>Sb + C„6 .!>6 



S.D 



- luih — h'nih - hnxh — 7„' dQ. 



(46) 



To evaluate the next three equations involving the fields, we make use of 
the fact that the expression for dU in equation (42) is a perfect differential. 
As a consequence there are relations between the partial derivatives, 
namely 

dj\n = dS . 
dff dS m 



dT m 
d8 n 

We note also that 



dEn 



dEn 
dff 



8Q 

d8„ 



(47) 



88 



TO \ 

n /8,<j 



= 4ir0 mn 



(48) 



where /3 is the so called "impermeability" matrix obtained fiom the dielectric 
matrix e nm by means of the equation 

0_ = { -^-. (49) 



where A is the determinant 



A = 



fll , 


Cl2 , 


€13 


Cl2 , 


£00 


eSJ 


C13 , 


€22 , 


C33 



(50) 



and A m,n the minor obtained by suppressing the mth row and nth column. 
The partial derivatives of the fields by the entropy can be written 

«M fo-ifr) «*-='-£") «--*■« <*> 

d(T /S.D dff /S.D 6 Off /S D 

where 9' S „ ,D is a pyroelectric constant measuring the increase in field required 
to produce a zero charge on the surface when a heat dQ is added to the 
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crystal. Since the voltage will be of opposite sign to the charge generated 
on the surface of the crystal in the absence of this counter voltage a nega- 
tive sign is given to q s n ' . 
Finally the last partial derivative 

de\ j l ae\ „ j 1 ae\ Jrt dQ 
da 



?) *i|?) e*-"?) <Q = § (52) 

T/S.D 6 Off / S,D U Off /S,D PC„ 



represents the ratio of the increase in temperature due to the added amount 
of heat dQ when the strains and electric displacements are held constant. 
It is therefore the inverse of the specific heat at constant volume and constant 
electric displacement per gram of material times the density p. Hence 
the ten equations of equation (44) can be written in the generalized forms 

T n = CnlSl + ^712^2 + CnlSa + Cn'iSi + C n 'sSb + C n't St 

—tini8i — h'^fa — h" n zh — y S n D dQ 

E m = —IllmSl — himSi — JUmSz ~ lt4 m Si — IlimSfi ~ llQmSf, + 47T/3ml5l 

+ 4*&k + 4tt/& 3 - qlf dQ (53) 

j/-. r\i S,Dp , S,D C I S,D C i . S,D c . S,D c i S,D C1 

dB = — G[7i Si +7 2 Si 4- 73 S 3 + 74 Si + 7b 05 + 7e o 6 ] 

.-.r S,D- 1 S,D~ 1 S,D ? i 1 *~ja 

— 0[?i 5i + g 2 5 2 + ? 3 63] + "7^ • 

11 = 1 to 6; m = 1 to 3 

If, as is usually the case with vibrating crystals the vibration occurs 
with no interchange of heat between adjacent elements dQ = and the 
ten equations reduce to the usual nine given by the general forms 

T n = CnlSl + C„2S2 + CnaS^ + C n iSi 4" CnfrSb + Cnt^e 
£m — —hmSl — ItemSi — ll3mS 3 ~ IlimSi ~ I'bmSi — h Gm Se 

+ tor&id! 4- 4r/£ 2 52 + 4x^6, . 



(54) 



In these equations the superscript a has been dropped since the ordinary 
constants are adiabatic. The tenth equation of (53) determines the increase 
in temperature caused by the strains and displacements in the absence of 
any flow of heat. 

If we introduce the expression of equations (53) into equation (42) the 
total energy of the crystal is per unit volume. 
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W = cfi'Sl + 2c?2%6 1 2 + 2&S& + 2^4^54 + 2cf5 <r 5x5- 6 + 2c??S 1 S fi 

+ C&SI + 2C2Y5 253 + 2^254 + 2c^'SoS b + 2^^ 2 5 6 
+^"33 ,ff 'S'3 + 2c 3 4 ff 5'354 + 2czi"S 3 S 6 + 2c36 SaSe 

+ f44 ff »S4 -f- 2f45 S4S5 + 2^46 O4O6 
-J-C66 $6 T 2^66 O6O6 

+W (55) 

-(tfuWi + 2Aw«i5 2 + 2J&8i5i + 2A; 4 «i5 4 + 2hi S SiS 6 + 2Aw8i5 fl ) 

— (2/?2l5=^l + 2&fcS| + 2//23^ 3 + 2//24«2^4 + 2/i««*S6 + 2/;265 2 5 6 ) 

-(2*SASi + 2/45352 + 2/45 3 S 3 + 2/?3 4 53>S4 + 2//3 5 5 3 5 6 + 2//3 6 535 6 ) 
-(2 T f lD 5 1< /g + 2 7 2 ,D ^g + 2yi D S»dQ 

+ ^"WS + 275 ,D 5 6 ^ + 2 7 f' D 5 e ^) 
+4ir[/3?iV 1 + 2j8ff «x*t + 2/JiYM, + /3fjT5l + SpSftrii + jSffts] 

-{2q\' D hdQ + 2 9 f'°MQ + 2gf D W0 + ^T • 

pC- u 

Equations (53) can be derived from this expression by employing the partial 
derivatives 

dU _. dU ,_ 3t/ , 

r » = aT B ; *--aC' *-•&© (56) 

The other form for writing the elastic, piezoelectric, pyroelectric and di- 
electric relations is to take the strains, displacements, and entropy as the 
fundamental variables and the stresses, fields and temperature increments 
as the dependent variables. If we develop them in terms of their partial 
derivatives as was done in (44), use the relations between the partial deriva- 
tives shown in equation (57). 

d8m = dSn . dS n _ da . db» = da_ (57) 

8T n dE m ' 80 dT n ' dQ dE m 

and substitute for the partial derivatives their equivalent elastic, piezo- 
electric, pyroelectric, temperature expansions, dielectric and specific heat 
constants, there are 10 equations of the form 
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s n = slfT! + s' n f r, + s*fT 3 + s E n ?T 4 + 5*i e r 8 + ,*f r, + d Q nl E i 

+ 5*Lj& + e |?£ 2 + ^£ 3 + rfrfe (58) 

47T 4lT 4tf 

dQ = B da = Q[cti T\ + a 2 T 2 + 0:3 ^3 + «4 Z4 + «5 r 6 + a 6 T 6 ] 

» = 1 to 6, m = 1 to 3 

The superscripts -E, 0, and T indicate respectively constant field, constant 
temperature and constant stress for the measurements of the respective 
constants. It will be noted that the elastic compliance and the piezo- 
electric constants d mn are for isothermal conditions. The a s constants are 
the temperature expansion constants measured at constant field, while the 
p T constants are the pyroelectric constants relating the ratio of 5 = D/4ir 
to increase in temperature dQ, measured at constant stress. Since there is 
constant stress, these constants take into account not only the "true" pyro- 
electric effect which is the ratio of 5 = Z>/47r to the temperature at constant 
volume, but also the so called "false" pyroelectric effect of the first kind 
which is the polarization caused by the temperature expansion of the crystal. 
This appears to be a misnomer. A better designation for the two effects 
is the pyroelectric effect at constant strain and the pyroelectric effect at 
constant stress. C% is the specific heat at constant pressure and constant 
field. 

If we substitute these equations into equation (42), the total free energy 
becomes 

w = E 2 <Er«r. + 2 £ S P«*tje* + 2 £ alrje 

m= 1 71—1 71=1 0«=1 71—1 

3 B (59) 

+ E i t^E E p + 2 ZpoE.de + & dQ. 
0=1 P =i 4tt o-=i fc» 

Equation (58) can then be obtained by partial derivatives of the sort 

_ HE <w_ _<& _ dU 

Sn " ar„ ; 8o " d£ P ; rf<r " e " d(<*e)' 

By tensor transformations the expression for U in (59) can be shown to 
be equal to the expression for U in (55). 
The adiabatic equations holding for a rapidly vibrating crystal can be 
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obtained by setting dQ equal to zero in the last of equations (58) and elim- 
inating dQ from the other nine equations. The resulting equations are 

S„ = s„iT\ + s n2 T2 + s% 3 T 3 + SniTi 

+ s*h 7"b + Sni T& + d n iEi + d n iEi + d ni Ei 

8 m = d lm Ti + d im T 2 + d 3m T 3 + d im T 4 (*) 

4tt 4tt 47r 

where the symbol a for adiabatic is understood and where the relations 
between the isothermal and adiabatic constants are given by 

E E r. E .T a T,a T,e .T .T f~. 

E,a _ E,Q _ CX m a„^ _ ja _ , _ <Xtn ${ tl € mn _ C mn _ pmpnU 

pC» ' ae ~ dlm ~pCT ; 4t " 4^ P Cf ' 

Hence the piezoelectric and dielectric constants are identical for isothermal 
and adiabatic conditions provided the crystal is not pyroelectric, but differ 
if the crystal is pyroelectric. The difference between the adiabatic and 
isothermal elastic compliances was discussed in section (1.4) and was shown 
to be small. Hence the equations in the form (60) are generally used in 
discussing piezoelectric crystals. 

Two other mixed forms are also used but a discussion of them will be 
delayed until a tensor notation for piezoelectric crystals has been discussed. 
This simplifies the writing of such equations. 

3. General Properties of Tensors 

The expressions for the piezoelectric relations discussed in section 2 can 
be considerably abbreviated by expressing them in tensor form. Further- 
more, the calculation of elastic constants for rotated crystals is considerably 
simplified by the geometrical transformation laws established for tensors. 
Hence it has seemed worthwhile to express the elastic, electric, and piezo- 
electric relations of a piezoelectric crystal in tensor form. It is the purpose 
of this section to discuss the general properties of tensors applicable to 
Cartesian coordinates. 

If we have two sets of rectangular axes (Ox, Oy, Oz) and (Ox', Oy', Oz) 
having the same origin, the coordinates of any point P with respect to the 
second set are given in terms of the first set by the equations 

x' = l\X + tthy + «i« 

y' = kx + nhy + «jz (61) 

%' = l%x + m%y + w»z. 
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The quantities (A , • • • , n 3 ) are the cosines of the angles between the various 
axes; thus li is the cosine of the angle between the axes Ox' , and Ox; n 3 the 
cosine of the angle between Oz' and Oz, and so on. By solving the equations 
(61) simultaneously, the coordinates x, y, z can be expressed in terms of 
x', y', z' by the equations. 

x = lix' + ky' + l 3 z' 

y = nti/x' + nny' + m 3 z' (62) 

z = n\x' + thy' + «3z'. 

We can shorten the writing of equations (61) and (62) considerably by 
changing the notation. Instead of x, y, z let us write X\ , x? , x 3 and in place 
of x', y', z' we write x\ , a* 2 , x 3 . We can now say that the coordinates with 
respect to the first system are Xi , where i may be 1, 2, 3 while those with 
respect of the second system are x, , where j = 1, 2 or 3. Then in (61) 
each coordinate x,- is expressed as the sum of three terms depending on the 
three x,- . Each Xi is associated with the cosine of the angle between the 
direction of Xi increasing and that of X, increasing. Let us denote this 
cosine by a,; . Then we have for all values of j, 

3 

x] = aijXi + a 2) x 2 + (hi x 3 = 2 0.7 *i- (63) 

»— i 

Conversely equation (62) can be written 

3 

7=1 

where the a,-, have the same value as in (63), for the same values of i a.ndj, 
since in both cases the cosine of the angle is between the values of Xi and Xj 
increasing. Such a set of three quantities involving a relation between two 
coordinate systems is called a tensor of the first rank or a vector. 

We note that each of the equations (63), (64) is really a set of three equa- 
tions. Where the suffix i or j appears on the left it is to be given in turn 
all the values 1, 2, 3 and the resulting equation is one of the set. In each 
such equation the right side is the sum of three terms obtained by giving j 
or i the values 1, 2, 3 in turn and adding. Whenever such a summation 
occurs a suffix is repeated in the expression for the general term as a,-,* ; - . 
We make it a regular convention that whenever a suffix is repeated it is 
to be given all possible values and that the terms are to be added for all. 
Then (63) can be written simply as 

Xj = aijXi 

the summation being automatically understood by the convention. 
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There are single quantities such as mass and distance, that are the same 
for all systems of coordinates. These are called tensors of the zero rank 
or scalars. 

Consider now two tensors of the first rank u ,- and z>* . Suppose that each 
component of one is to be multiplied by each component of the other, then 
we obtain a set of nine quantities expressed by u ,• 0* , where i and k are 
independently given all the values 1, 2, 3. The components of «,- v* with 
respect to the Xj set of axes are u, Vf , and 

u'j v'( = (anUi) (akfik) = an a*t una, . (65) 

The suffixes i and k are repeated on the right. Hence (65) represents nine 
equations, each with nine terms. Each term on the right is the product 
of two factors, one of the form a,;a^, depending only on the orientation of 
the axes, and the other of the form u&u , representing the products of the 
components referred to the original axes. In this way the various «,- zy can 
be obtained in terms of the original uflk . But products of vectors are not 
the only quantities satisfying the rule. In general a set of nine quantities 
Wit referred to a set of axes, and transformed to another set by the rule 

Wfi = anakiWik (66) 

is called a tensor of the second rank. 

Higher orders tensors can be formed by taking the products of more 
vectors. Thus a set of n quantities that transforms like the vector product 
Xi Xj • • • x p is called a tensor of rank n, where n is the number of factors. 

On the right hand side of (66) the i and k are dummy suffices; that is, 
they are given the numbers 1 to 3 and summed. It, therefore, makes no 
difference which we call i and which k so that 

Wjg = anakeWik = akjai(K'k( . (67) 

Hence Whf transforms by the same rule as w ik and hence is a tensor of the 
second rank. The importance of this is that if we have a set of quantities 



(68) 



Wn 


K'l2 


U<13 


Wn 


ft'22 


2£'23 


1031 


Wgt 


ft'33 



which we know to be a tensor of the second rank, the set of quantities 



tt'll 


"'21 


1031 


«'12 


W'22 


1031 


wis 


1023 


1033 



(69) 



is another tensor of the second rank. Hence the sum (wtk + «'*<) and the 
difference (wa — «'**) are also tensors of the second rank. The first of 
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these has the property that it is unaltered by interchanging i and k and 
therefore it is called a symmetrical tensor. The second has its components 
reversed in sign when i and k are interchanged. It is therefore an antisym- 
metrical tensor. Clearly in an antisymmetric tensor the leading diagonal 
components will all be zero, i.e., those with i = k will be zero. Now since 

wn, = \ (w ik + Wki) + \ (wik — ivki) (70) 

we can consider any tensor of the second rank as the sum of a symmetrical 
and an antisymmetrical tensor. Most tensors in the theory of elasticity 
are symmetrical tensors. 

The operation of putting two suffixes in a tensor equal and adding the 
terms is known as contraction of the tensor. It gives a tensor two ranks 
lower than the original one. If for instance we contract the tensor m Vk 
we obtain 

UiVt = uivi + U2V2 + U3V3 (71) 

which is the scalar product of u » and Vk and hence is a tensor of zero rank. 

We wish now to derive the formulae for tensor transformation to a new 
set of axes. For a tensor of the first rank (a vector) this has been given 
by equation (61). But the direction consines t\ to n 3 can be expressed in 
the form 



(72) 



dx' dxi 
dx dzi ' 


dx' dxi 
dy dX2 


dx' 

ft! = — 

dz 


3«i 
dx 3 


dx dxi ' 


dy' dx 2 
dy dxi 


dy' 
dz 


dx 2 
dx 3 


_ dz' _ dx' 3 
dx dxi ' 


dz' dx 3 
dy 0x2 


dz' 
dz 


__ dx 3 
dx 3 


Hence equation (61) can 


be expressed in the tensor form 






/ _ dx'j _ 

Xj — — Xi — UijXi . 

oXi 







(73) 

Similarly since a tensor of the second rank can be regarded as the product 
of two vectors, it can be transformed according to the equation 

/ / /dx'j \ (dx, \ dx, dx, , . 

which can also be expressed in the generalized form 

/ dx, dx/ • .„-. 

""-ai/d""*- (75) 
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In general the transformation equation of a tensor of the «th rank can be 
written 



, _ dXki dXk* 

A *>-*» a. Vyi dx i2 



fok'n Y 

^7 Ay,,,.../.. 



(76) 



4. Application of Tensor Notation to the Elastic, Piezoelectric 
and Dielectric Equations of a Crystal 

Let us consider the stress components of equation (7) 



1 yy A * 

7\„ T, 



from which equation (8) is derived 

I xy = I yx > I xz = I zx \ 1 yz = *■ zy 

and designate them in the manner shown by equation (77) to correspond 
with tensor notations 



T\\ 7*12 ^13 
Til T 22 T 23 
7*31 7 32 3- 83 



7*11 ^12 *13 
Ttf ^22 ^23 
jTi3 J 23 ^33 



(77) 



by virtue of the relations of (8). We wish to show now that the set of 9 
elements of the equation constitutes a tensor, and by virtue of the relations 
of (8) a symmetrical tensor. 

The transformation of the stress components to a new set of axes x', y', z' 
has been shown by Love to take the form 

T xx = l\ T xx + m\Tyy + n\T zz + 2twiT iy + 2CvhT X z + 2m x n{T yi 

T' XII = IxkT^ + nwnvTyy + WfoTn + (A«i + tim)T iy + (A»a + hni)T x , 

+ (wi»2 + taHh)T wn 

where A to « 3 are the direction cosines between the axes as specified by 
equation (61). Noting that from (72) 



(78) 



f _ dx i 
(l -dx~r 



dx 3 
OXa 



the first of these equations can be put in the form 

4 See "Theory of Elasticity," Love, Page 80. 
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,' ( dx '*\ rr, , dx[ dx[ , dx[ dx[ 

11 = I jr— I i n i - 7— ^ — In -r x— % — i 11 



dx\ 8x2 



dx[ dx[ , (dxT\ 

1 Tn+ \te,) 



dxi dxg 



dxi dx\ 



dxi dxi 



dx 2 dxi 



, ojl.1 0x1 0x1 qxi 

~r -r — ■= — 1 23 — r — •= — i */ 

0x2 dxz dXk oX( 



(79) 



, dxi dxi , dxi dx! , / dxi\ 

d.r 3 Mi d£ 3 0x2 \o%i/ 



while the last equation takes the form 



/ dxi dx 2 ~ , dxi dx 2 „ , dxi dx 2 - 

i 12 — t — t — 1 11 T r — -r— i 12 T ~ — -^ — i 13 
OXi OXi OXi 0X2 OXi 0X3 

, dx[ dx 2 rp , 5*i 3*2 d*i dx'2 rr, dx'i dx' 2 

T -5 — -r — i 21 i~ -t: — t — J- 22 ~t~ -r — - — i 23 — r — -r — i fc/ 

0X2 0x1 0x2 0x2 0x2 ox 3 dXh oxf 

, dxi dx'2 rr , dx[ dx'2 rr, , d*( 3*s -, 

~r -5 — ^ — J- 31 ~r t — t — i 32 ~r -5 — 5 — ■* 33 

0x3 dxi 6x3 0x2 0x3 0x3 

The general expression for any component then is 

/ _ dx'i dx'j 
dx k dx f 



(80) 



(81) 



which is the transformation equation of a tensor of the second rank. Hence 
the stress components satisfy the conditions for a second rank tensor. 
The strain components 



s xt 


^xy 


s z . 


Oyx 


Syy 


Of/t 


s zx 


lJ zy 


s„ 



do not however satisfy the conditions for a second rank tensor. This is 
shown by the transformation of strain components to a new set of axes, 
which have been shown by Love to satisfy the equations 

S xx = l\ S xx + miSyy + n\S z , + ImiSxy + liniS zx + m\n\S V z 



Sxy = 2tit 2 S xx + 2m\tn2Syy + 2miw 2 'Sm + {l\m 2 + hmijSxy 

+ (Aws + niQSn + (wiM 2 + ffhn^Sx, . 
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If, however, we take the strain components as 



«? - 9 J d * ■ s 

OH — Oix — — , ^ 


22 ~ 6w ~ a? ' M " " ~ dz 


e c - Sx « - 1 ( dr l 4- d t\ • c _ C - S " 


1 ^ 4. ^ • C - 9 - S " - l ( d * 4- d A 

- 2 \Fx + ax) • 523 - 5a2 - t " 2 U *; 


the nine components 




Si} S\2 Ol3 






521 022 023 






•S31 6^32 033 





(83) 



(83) 



will form a tensor of the second rank, as can be shown by the transformation 
equations of (82). 
The generalized Hooke's law given by equation (22) becomes 



Ta = CijkiSki 



(84) 



dikl is a fourth rank tensor. The right hand side of the equation being 
the product of a fourth rank tensor by a second rank tensor is a sixth rank 
tensor, but since it has been contracted twice by having k and I in both 
terms the resultant of the right hand side is a second rank tensor. Since 
dm is a tensor of the fourth rank it will, in general, have 81 terms, but on 
account of the symmetry of the T ti and S k( tensors, there are many equiva- 
lences between the resulting elastic constants. These equivalences can be 
determined by expanding the terms of (84) and comparing with the equiva- 
lent expressions of (22). For example 



T\\ = cmio , i2 + cimo^ + C1113S13 

~f" CftftSn + Cji22«S'22 ~T" C1123023 
+ CmiSzl + Cu32/Ss2 + C1133S33 . 



(85) 



Comparing this equation with the first of (22) noting that Su — S21 — 



-y , etc., we have 



C\\w = Cn ; Ciiu = Cu2i = Cie ; cnas — Cis ; cnu — disi — £i6 \ 

Cll22 = £l2 J CUM = C U32 = c u • 



(86) 
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In a similar manner it can be shown that the elastic constants of (22) 
correspond to the tensor elastic constants cakt according to the relations 

Cll — C\n\ ', C12 = C1122 = £2211 ! Cl3 — C;i33 = £3311 ', Cl4 = ^1123 = ^1132 = 
C2311 = C3211 ' p c\b = Cnia = c 1131 = ^1311 = Qui ', ClS = Clll2 = Cll21 = Cl211 = 

^2111 ; C22 = C2222 j c<a ~ C2233 = ^ 3322 » C2i = C2223 = C2232 = f2322 = C3222 i 

C26 = ^2213 = ^2231 = ^1322 — C3122 ; ^26 = C22 l 2 = C2221 = C 1222 = f2122 > f33 = 
^3333 I £34 = c 3323 — C3332 = C2333 = ^3233 J £35 = 3313 = ^3331 = ^1333 = C3133 J 

(87) 

^30 = c 3312 — £3321 — ^1233 — C 2133 J £44 — ^2323 — C 2332 — ^3223 — £"3232 5 ^46 — 
i"2313 = C2331 = ^S 213 = C3231 = 1323 = 133Z = C3132 = C3123 » Ci6 = C2312 = 
<"2321 = £3212 = ^3221 = ^1223 = ^1232 = ^2123 = ^2132 ', ^66 = C1313 = £1331 = 
C3113 = C313I J ^86 = C1312 = fl321 = ^3112 = C3121 = C1213 = <"l231 = ^2113 = 
C2131 J CG6 = Cl212 = Cl221 = f2112 = C2121 ' 

Hence there are only 21 independent constants of the 81 a,hg constants 
which are determined from the ordinarily elastic constants c.-y by replacing 

1 by 11; 2 by 22; 3 by 33; 4 by 23; 5 by 13; 6 by 12 (88) 

and taking all possible permutations of these constants by interchanging 
them in pairs. 
The inverse elastic equations (26) can be written in the simplified form 

Sij= s iik( T k( . (89. 

By expanding these equations and comparing with equations (26) we can 
establish the relationships 

5 " _ 
Sn = Sun ; Su = S1122 — *aii 5 5i3 — ^1133 — 53311 ; ~r — ^1123 — 51132 — $2311 — 

5l6 Sis _ 
53211 J -T = ^1113 — ^1131 — 51311 — 53111 ; ~z 5ni2 — 5ll2i — 5l2ll — 52111 ', 

_ 524 _ 526 

522 = 52222 ', 523 — 52233 — 53322 ; — — 52223 _ 52232 — 52322 — 53222 ; ~T — 

(90 A) 

Sn _ _ __ _ 

52213 — 52231 — 5l322 — 53122 ', — — 52212 — 52221 — 5i222 — 52122 ', S33 — S3333 ; 

S M _ „ Ssl — c — * _ 536 _ 

— — 53323 — 53332 — 52333 — 53233 ; ~Z ~ 53313 — $3331 — 51333 — S3133 , — 
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__ $44 _ _ _ S45 

$3312 — S3321 — ^1233 — $2133 5 -T — $2323 — fan ~ $3223 — $3232 ', ~r = fal3 = 

•^2331 — $3213 — $3231 ~ ^1323 — $1332 — $3123 = $3132 5 — = $2312 = ^2321 = 

(90 B) 

$55 
$3212 — -S3221 — $1223 — $1232 — $2123 — $2132 ', ~T — $1313 — $1331 — $3113 = 

$66 _ 
$3131 J -7 $1312 — $1321 — $3112 — $3121 = $1213 — $1231 = $2113 = $2131 ', 

$66 _ 

-7 $1212 — $1221 — $2112 — $2121 • 

Here again the Sijtf elastic constants are determined from the ordinary 
elastic constants stf by replacing 

1 by 11, 2 by 22, 3 by 33, 4 by 23, 5 by 13, 6 by 12. 

However for any number 4, 5, or 6 the elastic compliance s {i has to be di- 
vided by two to equal the corresponding suit compliance, and if 4, 5 or 
6 occurs twice, the divisor has to be 4. 

The isothermal elastic compliance of equations (39) can be expressed 
in tensor form 

Sn= s%ktT k t+aijde (91) 

where as before a*, is a tensor of the second rank having the relations to 
the ordinary coefficients of expansion 

ot\ = «n , oil — oiii , at = ass J ■=- — 0:23 i 

y = «13 J y = ai2. 

The heat temperature equation of (35) is written in the simple form 

dQ = + a h t T k t G + pC p dB. (92) 

By eliminating dQ from (92) and substituting in (91) the adiabatic constants 
are given in the simple form 

Sijkl — s ijlc e — — - — . (93) 

pC p 

The combination elastic and piezoelectric equations (60) can be written 
in the tensor form 

S tj = s E i]k( T k( + d mij E m ; 8 n = e -^E m + d nk (T k (. (94) 



■ 
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Here d m a is a tensor of third rank and Cm„ one of second rank. The </„,,-, 
constants are related to the eighteen ordinary constants da by the equations 



d\\ = dm ', d\i — dy 



di3 = dia3 ; — = di2 3 = dm ; -= dna — dm 



du 

2 



-r^ = dm = dm ', ^21 = dm ', dyi = d-m ; d^z — d<x3 ; — d>m — <mi 



2 



-^ = dm = dl3\ ! -=" = ^212 = ^221 J ^31 = ^311 \ ^32 — ^322 i ^33 — ^333 



(95) 



" = ^323 — ^332 J -=" = ^313 — ^331 \~Z '' r - i 



(1st 

2 



U ' 36 A - i 

— i':wj. - <<;•; 

m 



The tensor equations (94) give a simple method of expressing the piezo- 
electric equations in an alternate form which is useful for some purposes. 
This involves relating the stress, strain, and displacement, rather than the 
applied field strength as in (94). To do this let us multiply through the 
right hand equation of (94) by the tensor 47r/3 mn , obtaining 

AlT&lnbn = tZn&nEn + tod nh l fc n T k t (96) 

where /3 m „ is a tensor of the "free" dielectric impermeability obtained from 
the determinant. 



Mmn 



whe- e A* is the determinant 



= (-D 



(m+n) 



n) &mn 



A' = 



T 

til 



111 



T T 

612 *22 



T T 

«13 *23 



T 

CIS 



T 
621 



<33 



(97) 



(98) 



and Amn the minor obtained from this by suppressing the f»th row and nth 
column. If we take the product <= mn /3 mn for the three values of m, we have 
as multipliers of Ei , E 2 , E 3 , respectively 



en /3n + €12 pit -T ei3 Pis — 1 

«21 021 + €22 #22 + «23 #23 = 1 
€31 P31 + €32 P82 T C33 P33 — 1. 



(99) 



But by virtue of equations (97) and (98) it is obvious that the value of 
each term of (99) is unity. Hence we have 



£,» = 4r/3 ra „ 8 n - (4tt d nk ( m „) T k e . 



(100) 



PIEZOELECTRIC CRYSTALS IN TENSOR FORM 



109 



Since the dummy index n is summed for the values 1, 2, and 3, we can set 
the value of the terms in brackets equal to 

gmkl = 47T dnkt Pmn = Wfati fill + fatf fljfl + d Sk ( ftj,] (101) 

and equation (100) becomes 



E m = 47r fi mn 8 n — g mk g Tut . 



(102) 



Substituting this equation in the first equations of (94) we have 

S ii = s? ik (T k i + g; ii 8 n (103) 

where 

SijkC = Sijkt — d mij g mk e = s* ik i — 47r[/3mn d nk ( dmij]. 

By substituting in the various values of i,j, k and I corresponding to the 21 
elastic constants, the difference between the constant displacement and 
constant potential elastic constants can be calculated. If equations (102) 
and (103) are expressed in terms of the S u -- ■, S 6 strains and Zi,-*-, T t 
stresses, the g„ ,-, constants are related to the g i} - constants as are the corre- 
sponding dij constants to the d nij constants of equation (95). 

Another variation of the piezoelectric equations which is sometimes em- 
ployed is one for which the stresses are expressed in terms Of the strains 
and field strength. This form can be derived directly from equations (94) 
by multiplying both sides of the first equation by the tensor c E ijk t for the 
elastic constants, where these are defined in terms of the corresponding 
Si,- elastic compliances by the equation 

4 = (-l) C,+y) A:-?/A j£ (104) 

where A is the determinant 

S B B B E B 

$11 $12 $13 $14 $16 $16 

B B E E B B 

$12 $22 $23 $24 $25 $26 

E B B BE E 

$13 $23 $33 S34 $35 $36 

B B B E B B 

$14 $24 $34 $44 $45 $48 

B B E E E B 

$J6 $26 $36 $45 $ 65 $ 66 

M E E E E B 

516 $26 $36 $46 $56 $66 

and A',- in the minor obtained by suppressing the *th row and^'th column. 
Carrying out the tensor multiplication we have 



A'* = 



Cij k (. Sn = c E ijk c s S ij k t T k e + d mij c B i ik t E„ . 



(105) 
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As before we find that the tensor product of c B ijk t s%ht is unity for all values 
of k and I. Hence equation (105) can be written in the form 

T k t=c E iik tSii-e m ucE m (106) 

where e mk ( is the sum 

e mk t = dmij c E iik t ( 107 ) 

EumxE ed for all values of the dummy indices i and j. If we substitute the 
equation (106) in the last equation of (94) we find 

5 n = ^E m + e nij S u (108) 

where c™„ the clamped dielectric constant is related to the free dielectric 
constant el„ by the equation 

e s mn = Zn-±ir[dn k te mk t}. (109) 

Expressed in two index piezoelectric constants involving the strains Sn m • -Su 
and stresses 7V • • T u the relation between the two and three index piezo- 
electric constants is given by the equation 

en = em ] «ie = e m J «n = £133 ; «M — g i23 = ftn ; en = e m = e m 
eie = em = em ; c\\ = e m ; e 2 2 = em ; «23 = em \ eu = £223 " g 232 

C2B = 0213 = ^281 ! g26 = em ~ em ' eU = em ' e3i = em ' ^ = 6m 

C34 = ^323 = 6332 ; e 3 & = em = e 33 i ; e 3 & = e 3 u = £321 • 

Finally, the fourth form for expressing the piezoelectric relation is the 
one given by equation (53). Expressed in tensor form, these equations 
become 

T k t = fasti - J', k C 5„ ; E m = 4qf&. 5„ - kmii'Su (111) 

In this equation the three index piezoelectric constants of equation (111) are 
related to the two index constants of equation (53) as the e constants of 
(110). These equations can also be derived directly from (106) and (108) 
by eliminating E m from the two equations. This substitution yields the 
additional relations 

h, lk t = ^Tre mk tf3mn ; C°j k t = Cijkt + 6mkl ^mij = Cijkl 

+ 4tt emicC e, t ij /3 m „ 

where 

&„ = (-l) (m+n) A; S „/A <S 



(110) 



(112) 
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in which 



8 

en 


.s 

eia 


S 
613 


a 


s 

622 


S 
623 


a 
613 


5 
623 


633 



A'* = 



The four forms of the piezoelectric equations, and the relation between 
them are given in Table I. 

Table I 

Four Forms of tiie Elastic, Dielectric, and Piezo Electric Equations 

and their interrelations 



Form 


Elastic Relation 


Electric Relation 


1 


Sij = Sijk(Tkl + dmijEm 


5„ 


T 


2 


Sij = SijktTkf + gnij&n 


£m 


= 47Tpl„S n - g m klTkt 


3 


T H = CiilkSij — emkCEm 


5n 


-£s*. + ..,*„ 


4 


Tkt = CijkfSij — llnkfSn 


E m 


= 4irp S mn 8 n - KijSij 



Form 


Relation Between 
Elastic Conaliin.^ 


Relation Between 
Piezoelectric Constants 


Relation Between 

Dielectric Constants 


1 


sD iikt =S lkt- d ^M 


gmkl= 4,1-/3^ d nk( 


f 3 L = (- i y m+n)A 'J n / A ' T 


2 


c*. = (-l) ii+i) A'.*/A- E 


e mk l = *miicf fk( 


e ln = «L ~ *K**ft*l) 


3 


'me^'me+'-rf-u 


Kkl = W s mn e mkf 


OS _oT I gnkf.l'mkl 
P mn P mn ' 4„. 


4 


c o_ = (-l)^ + >) A ^/A sD 


Kkt = gniiC D iikt 


^ n= (_l)<" + »> A £/ A eS 



5. Effect of Symmetry and Orientation on the Dielectric Piezo- 
electric and Elastic Constants of Crystals 

All crystals can be divided into 32 classes depending on the type of sym- 
metry. These groups can be divided into seven general classifications 
depending on how the axes are related and furthermore all 32 classes can 
be built out of symmetries based on twofold (binary) axes, threefold (trig- 
onal) axes, fourfold axes of symmetry, sixfold axes of symmetry, planes of 
reflection symmetry and combinations of axis reflection symmetry besides 
a simple symmetry through the center. Each of these types of symmetry 
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result in a reduction of the number of dielectric, piezoelectric, and elastic 
constants. 

Since the tensor equation is easily transformed to a new set of axes by 
the transformaion equations (76) this form is particularly advantageous 
for determining the reduction in elastic, piezoelectric and dielectric con- 
stants. For example consider the second rank tensors tU and atf for the 
dielectric constant and the expansion coefficients. Ordinarily for the most 
general symmetry each tensor, since it is symmetrical, requires six inde- 
pendent coefficients. Suppose however that the X axis is an axis of twofold 
or binary symmetry, i.e., the properties along the positive Z axis are the 
same as those along the negative Z axis. If we rotate the axes 180° about 
the X axis so that + Z is changed into - Z, the direction cosines are 

. 9x[ . dx[ n dxi _ „ 

ft ■ — = 1 i Wi = — = ; n t = — =V 

1 dxi dx 2 dx 3 

a dx 2 n dx * 1 . « — dXi — n (WW 

dxi dx 2 ox a 

. dx 'i n dx * - n. « - dx ' 3 - -1 

The tensor transformation equations for a second rank tensor are 

, _ dx'i dx'i au x 

dXk dxi 
Applying (113) to (114) summing for all values of k and / for each value of 
i, and j we have the six components 

€11 = «11 5 «12 = — C 12 5 €13 = — €13 ', €22 = €22 5 € 2 3 = €23 \ €33 = 633 • (H5) 

Since a crystal having the A' axis a binary axis of symmetry must have the 
same constants for a +Z direction as for a — Z direction, this condition 
can only be satisfied by 

€12 = €13 = 0. (116) 

The same condition is true for the expansion coefficients since they form a 
second rank tensor and hence 

an = on = 0. (117) 

In a third rank tensor such as di ik , em, , gnu , htjh , we similarly find that 
of the eighteen independent constants 

hm = hi ; //lw = hi ; //2u = ha. ; A222 — fa ; hm = Am ; 

(llo) 

hm = ha ; h 3 n = /'31 ; /'322 = I132 ; /'323 = h 3 i ; hm = A33 • 
are all zero. The same terms in the d iih ,«»,*, giih tensors are also zero. 



PIEZOELECTRIC CR YSTA LS IN TENSOR FORM 1 1 3 

In a fourth rank tensor such as c ,-,*/, s iik t, applying the tensor trans- 
formation equation 

dXi dxj dxt dxl 
and the condition (113) we similarly find 

CIS = Ci6 = C 2 B = C 26 = C 36 = C 3 « = C45 = C« = 0. (120) 

If the binary axis had been the Y axis the corresponding missing terms 
can be obtained by cyclically rotating the tensor indices. The missing 
terms are for the second, third and fourth rank tensors, transformed to 
two index symbols, 

623 , €12 ; An , A12 , hm , hi , fiu , A 26 , h 3 i , A 32 , A33 , A35 ; 

(121) 

Cli , Ci« , C24 , C26 , £34 , C36 , C45 , C55 • 

Similarly if the Z axis is the binary axis, the missing constants are 

cis , €12 ; hu , h\2 , ht» , h\6 , hn , hyi , Ihz , /^e , A34 , A«6 ; 

(122) 

Cl4 , Cib , C24 , C25 , Cm , C3B , C46 , C B 6 • 

Hence a crystal of the orthorhombic bisphenoidal class or class 6, which 
has three binary axes, the X, Y and Z directions, will have the remaining 
terms, 

en , €22 , €33 ; hu , hit, , has ; c\\ , Cn , C13 , c 22 , C23 , C33 , Cu , c 65 , c 6t (123) 

with similar terms for other tensors of the same rank. Rochelle salt is a 
crystal of this class. 

If Z is a threefold axis of symmetry, the direction cosines for a set of 
axes rotated 120° clockwise about Z are, 

/ d#i c dx i £ti d x 'i 

l\ = r— = - .5 ; mi = — = - .866 ; n x = ~ = 
0x1 dx 2 dx 3 

(2 = ^ = .866; m 2 =^=-.5; n 2 = ^ = (124) 

0X1 dxz dx 3 

h = — = ; m s = — = ; » s = _£ = 1. 

o»i d# 2 ds 3 

Applying these relations to equations (114) for a second rank tensor, we 
find for the components 

en = .25 €n+ .433 en 4- .75e22 ; €12 = — .433 e u + .25ei 2 + .433 e 22 

di = —.5613 —.866623; €22 = .75en — .433ei2 + .25e2» (125) 

€23 = .866611 — -562J ; 631 = 6|3 . 
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For the third and fifth equations, since we must have ei 3 = €13 ; e 2 3 = fta 
in order to satisfy the symmetry relation, the equations can only be satis- 
fied if 

Q| = €23 = 0. (126) 

Similarly solving the first three equations simultaneously, we find 

€12 = 0; «U= C22- ( 127 ) 

Hence the remaining constants are 

en = C22 ; f33 • (128) 

Similarly for third and fourth rank tensors, for a crystal having Z a trigonal 
axis, the remaining terms are 

#11 , ha = — hi , hz = 0; #u , h , ^ie = — ha 

} l21 = -i h2 , h<n , #23 = 0, h u = h ; #25 = -*u , hi = -hi (129) 

//31 ; #32 = #31 ; #33 ; #34 = 0; #35 = 0; #36 = 

di ; C\2 ; C13 ; Cu ', ci = — ca] tie = 

C12 ; £22 = c u ; c 2 3 = C13 ; C24 = — cu ; C25 ; Cae = 

c\z ; to = C13 ; C33 ; C34 = 0; c 35 = 0; c a6 = 

C14 ; C24 = — cu ; C34 = 0; cu ; C45 = 0; C46 = cib 

Ci6 = — C2B ; C26 ; C35 = 0; C45 = 0; c 6 6 = cu ; c 00 = cm 

Cie = 0; C26 = 0; c 3 6 = 0; c i6 = c 2o ; c b <$ = Cm ; c 66 — \ (cu-Cjs). 

If the Z axis is a trigonal axis and the X a binary axis, as it is in quartz, 
the resulting constants are obtained by combining the conditions (116), 
(118), (120) with conditions (128), (129), (130) respectively. The resulting 
second, third and fourth rank tensors have the following terms 

«ii ; €12 = 0; €13 = 

€12 = 0; €22 = en ; «23 = (131) 

€13 = 0; €23 = 0; €33 

hi ; hi = -hi ; #13 = 0; hu ; ho = 0; he = 

#21 = 0; ha = 0; hi = 0; # 2 4 = 0; hb = ~hi ;ho = ~hi (132) 

#31 = 0; #32 = 0; #33 = 0; #34 = 0; #35 = 0; # 36 •= 
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(133) 



Cu ', Cn ', C13 ; Cu ; C15 = 0; Ci8 = 

Cn ', c-ii = C\\ ; Ciz = c\% ; c»a = — Cu ', c%i = 0; c^ = 

c\3 ; C23 = Cn ; C33 ; C34 = 0; C35 = 0; C36 = 

Cm ; C14 — — C14 ; r 3 4 = 0; Cu ; C45 = 0; C46 = 

^15 — 0; r 2 5 — 0; C35 = 0; C45 = 0; C55 = C44 ; C56 = Cu 

fig = 0; c 2 6 = 0; c 38 = 0; r 38 = 0; c 66 = c i4 ; c 66 = £ (cii-tfia). 

5.1 Second Rank Tensors for Crystal Classes 

The symmetry relations have been calculated for all classes of crystals. 
For a second rank tensor such as e,y, the following forms are required 

Triclinic Classes 1 and 2 e u , t 12 , e 13 

Cl2 , «22 , «23 
«13 , «23 , €33 

«ii , , eu 

, 622 , 
«13 , , €33 
€11,0 ,0 

, coo , 1 ,-; ! 

,0 , €33 
«i,0 ,0 
, «u , 

,0 , e 33 
cii,0 ,0 
, en , 
0,0, €,, 
5.2 Third Rank Tensors of the Piezoelectric Type for the Crystal Classes 



Monoclinic sphenoidal, Y a binary axis, Class 3 
Monoclinic domatic, Y a plane of symmetry, Class 4 
Monoclinic prismatic, Center of symmetry, Class 5 

Orthorhombic 
Classes 6, 7, 8 



Tetragonal, Trigonal 
Hexagonal 
Classes 9 to 27 

Cubic 

Classes 28 to 32 



Triclinic Assymetric (Class 1) No 
Symmetry 



hn , /*i2 , /J13 , hu , //15 , Aie 
hn , h 2 2 , #23 , hu , hn , hit 

#31 , #32 , #33 , hu , #35 , #33 
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Monoclinic Sphenoidal (Class 3) Y is 
binary axis 



Monoclinic domatic (Class 4) Y plane 
is plane of symmetry 



Triclinic pinacoidal, (center of symmetry) h = (Class 2) 

,0 ,0 , hu , , h u 
hi , hn , h-a , , ^26 , 
,0 ,0 ,/?84,0 ,h K 
hu , hu , hu ,0 , ^i6 , 

,0 ,0 , //24 , ,h» 

hi , hsz , hu ,0 , h 36 ,0 

Monoclinic prismatic (center of symmetry) h = (Class 5) 

Orthorhombic bisphenoidal (Class 6) ,0 ,0 , hu , ,0 
X, Y, Z binary axes Q Q Q Q )h ^ 

,0 ,0 ,0 ,0 ,h 36 

,0 ,0 ,0 ,h n ,0 

,0 ,0 , // 2 4 , , 

//31 , //32 , />33 , ,0 ,0 

Orthorhombic bipyramidal (center of symmetry) h = (Class 8) 

, ,0 , hu,h lb ,0 



Orthorhombic pyramidal (Class 7) Z 
binary, X, Y, planes of symmetry 



Tetragonal bisphenoidal (Class 9) 
Z is quaternary alternating 



Tetragonal pyramidal (Class 10) Z 
is quaternary 



Tetragonal scalenohedral (Class 11) Z 
quaternary, X and Y binary 



Tetragonal trapezohedral (Class 12) 
Z quaternary, X and Y binary 



, 0,0, -//is , hi , 
hi , -hi , , , , foe 
,0 ,0 , hu , hi , 
,0 ,0 , h 6 , -h u , 

hi,hi,h*,o , o ,o 

,0 ,0 ,h u ,0 ,0 
,0 ,0 ,0 , hu , 
,0 ,0 ,0 ,0 ,h 6 

,0 ,0 ,hu, ,0 
0,0,0,0, -hu , 
,0 ,0 ,0 , ,0 



(135) 
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Tetragonal bipyramidal (center of symmtery) A = (Class 13) 

,0 ,0 ,0 ,// 15 ,0 
,0 ,0 ,A 16 ,0 ,0 



Ditetragonal pyramidal (Class 14) Z 
quaternary, X and Y planes of 
symmetry 



Ihi , h n , h» , ,0 ,0 

Ditetragonal bipyramidal (center of symmetry) A = (Class 15) 

Trigonal pyramidal (Class 
16) Z trigonal axis 



Trigonal trapezohedral (Class 
18), Z trigonal, X binary 



Trigonal bipyramidal (Class 
19), X trigonal, plane of 
symmetry 

Ditrigonal pyramidal (Class 
20) Z trigonal, I' plane of 
symmetry 



Ditrigonal bipyramidal (Class 
22) Z trigonal, Z plane of sym- 
metry and Y plane of symmetry 

Hexagonal pyramidal (Class 23) 
Z hexagonal 



Hexagonal trapezohedral (Class 
24) Z hexagonal, X binary 





An, 


-Aii , 


> Am , 


A15 


— A22 


-Am, 


An , 


, A 15 . 


-Am 


-An 


h\ , 


A31 , A 33 


,0 , 








7) center of symmetry, A 


= 




An, 


-A u , 


. Au , 


, 





, 


,0 


, 


-Aii , 


-An 


o , 


,0 


,0 , 


, 





ha , 


-An,0 


,0 , 





-A 22 


-km, 


A22 , 


,0 , 





—An 


, 


,0 


, , 








o , 


,0 


, 


AlB 


— A22 


— Am , 


A ?2 , 


, A15 , 








Ihi , 


A31 , A33 


, 








21) center of symmetry, 


A = 


) 




An , 


-An , 


,o , 










, 


,0 


, 





-Aii 




o , 


,0 


,0 , 










o , 


,0 


, A« , 


Au 







o , 


,0 


> Am , 


-A M 







An, 


An , A33 


,0 , 










, 


,0 


, Am , 


, 







o , 


,0 


,0 , 


-Am, 







o , 


,0 


,0 , 
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Hexagonal bipyramidal (Class 25) center of symmetry, h = 



Dihexagonal pyramidal (Class 26) X 
hexagonal Y plane of symmetry 



,0 ,0 ,0 ,//i5,0 
,0 ,0 ,// 1B ,0 ,0 
Aa,A»,Ais,0 ,0 ,0 



Dihexagonal bipyramidal (Class 27) center of symmetry, h = 



Cubic tetrahedral-pentagonal-dedo- 
cahedral (Class 28) X, Y, Z binary 



,0 ,0 ,/ri4,0 ,0 
,0 ,0 ,0 , hu , 
,0 ,0 ,0 ,0 ,hu 



Cubic pentagonal-icositetetrahedral (Class 29) h = 

Cubic, dyakisdodecahedral (Class 30) center of symmetry, h = 



Cubic, hexakistetrahedral (Class 31) 
X, Y, Z quaternary alternating 



,0 ,0 ,//i 4 ,0 ,0 
,0 ,0 ,0 : //i4,0 
,0 ,0 ,0 ,0 ,hu 



Cubic, hexakis-oclahedral (Class 32) center of symmetry, h = 

This third rank tensor has been expressed in terms of two index symbols 
rather than the three index tensor symbols, since the two index symbols 
are commonly used in expressing the piezoelectric effect. The relations 
for the // and e constants are 



hu , hi 6 , hit are equivalent to hm , haa , hiu 



(136) 



in three index symbols, whereas for the </,•/ and ga constants we have the 
relations 



Ou a* a^ equiva ient to dm , d iU , d iu (137) 

2 ' 2 ' 2 4 

Hence the d t relations for classes 16, 18, 19, and 22 will be somewhat dif- 
ferent than the // symbols given above. These classes will be 
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Class 16 



Class 18 



Class 19 



Class 22 



du — dn d u dn — 2d» 2 
— di% dn d\z —du —2dn 
d 3l d 3l d 33 

du -d n d u 

-du -2d u 



dn -d n -2d 22 
-d 22 da -2d u 


dn ~dn 
-2d n 





(138) 



5.3 Fourth Rank Tensors of the Elastic Type for the Crystal Classes 



Triclinic System 
(Classes 1 and 2) 21 
moduli 



Monoclinic System 
(Classes 3, 4 and 5) 12 
moduli 
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Rhombic System 

(Classes 6, 7 and 8) 
9 moduli 



Tetragonal system, Z 
a fourfold axis (Classes 
9, 10, 13) 7 moduli 



Tetragonal system, Z a 
fourfold axis. X a two- 
fold axis (Classes 11, 
12, 14, 15) 6 moduli 



Trigonal system, Z a 
twofold axis, (Classes 
16, 17) 7 moduli 
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The s tensor is 
entirely analo- 
gous 



The s tensor is 
entirely analo- 
gous 



The s tensor is 
entirely analo- 
gous 







C-ib Cu 



Cll — C12 



The s tensor is 
analogous ex- 
cept that $4« = 

2">25 , $66 = 2Si4 , 
5 6 8 = 2(Sii — 5i 2 ) 
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Trigonal system, Z a 
trigonal axis, X a 
binary axis (Classes 
18, 20, 21) 6 moduli 



Hexagonal system, Z a 
sixfold axis, X a two- 
fold axis (Classes 19, 
22, 23, 24, 25, 26, 27) 
5 moduli 



Cubic system (Classes 
28, 29, 30, 31, 32) 3 
moduli 



Isotropic bodies, 2 

moduli 
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The s tensor is 
analogous ex- 
cept that 5 5 « = 

2Si4 , Sflfl = 

2 (5 U — su) 
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The s tensor is 
analogous ex- 
cept Sit = 
2(s n — s u ) 



The s tensor is 
entirely analo- 
gous 



The 5 tensor 
analogous ex- 
cept last three 
diagonal terms 
are 2(sn — Su) 






122 BELL SYSTEM TECHNICAL JOURNAL 

5.4 Piezoelectric Equations for Rotated Axes 

Another application of the tensor equations for rotated axes is in deter- 
mining the piezoelectric equations of crystals whose length, width, and thick- 
ness do not coincide with the crystallographic axes of the crystal. Such 
oriented cuts are useful for they sometimes give properties that cannot be 
obtained with crystals lying along the crystallographic axes. Such proper- 
ties may be higher electromechanical coupling, freedom from coupling to 
undesired modes of motion, or low temperature coefficients of frequency. 
Hence in order to obtain the performance of such crystals it is necessary to 
be able to express the piezoelectric equations in a form suitable for these 
orientations. In fact in first measuring the properties of these crystals a 
series of oriented cuts is commonly used since by employing such cuts the 
resulting frequencies, and impedances can be used to calculate all the pri- 
mary constants of the crystal. 

The piezoelectric equations (111) are 

Tk[ = Cijk(Sij — hnkfin ', E m = ^Pmn^n ~ llmijSij . (HI) 

The first equation is a tensor of the second rank, while the second equation is 
a tensor of the first rank. If we wish to transform these equations to another 
set of axes x', y', z' , we can employ the tensor transformation equations 

/ = dXk dX( T = dx^dX( 
H dx k dx t H dxkdx e 

• [ciiklSn + 2cnkCSu + 2cizhtSa + &h(Sn 

+ 24*5. + AuSmA -PP- fat*i + *bA + fctflll (140) 

dx k oxf 

E' m = 47T ^ [A fix + At it + A* J - ^ 

dx m ox m 

• [limuSu + 2h m \iSn + 2// m i3.S'i3 + hm-nS^ + 2// m 23'S'23 + h m 3iSaa\. 

These equations express the new stresses and fields in terms of the old strains 
and displacements. To complete the transformation we need to express 
all quantities in terms of the new axes. For this purpose we employ the 
tensor equations 

dxtdXi , .**"*' (uU 

; dXi dxj dx n 

dXi 
where tH f are the direction cosines between the old and new axes. It is 

OXi 

- i - 

obvious that ~ = ^-A and the relations can be written 
dXi ox i 
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_ dx[ _ dxi _ d*i _ dxi 

Cl = ait " 53 ; ' 2 " a*; ; ' 3 ~ d* 3 

d#2 5^2 d#2 

Wi = t - / ; f»j = t - / ; W3 = r~7 (142) 

0.T1 0.^2 aa"3 



d*3 dz 3 d* 3 

7 

3 



wi = . / ; «2 — • , / ; «3 = . 

0x1 ox 2 dx 



Hence substituting equations (141) in equations (140) the transformation 
equations between the new and old axes become 

/ o dx k dx f dXi dx, , dx' k dx' t dx n , 

dXk ox 1 dxi oxj dXk ox/ dx n 

, , (143) 

„/ . oS ox m dx n / dx m dXi dXj , 

d# m OXn OX m OXi OXj 

These equations then provide means for determining the transformation of 
constants from one set of axes to another. 

As an example let us consider the case of an ADP crystal, vibrating longi- 
tudinally with its length along the Xi axis, its width along the x 2 axis and 
its thickness along the x 3 axis, which is also the x 3 axis, and determine the 
elastic, piezoelectric and dielectric constants that apply for this cut when 
x[ is 6 = 45° from x\ . Under these conditions 

. dx[ dxi . dx'i dxi . 

t x = _=—,= cos 0; I, =-—=_, = -sm 0; 
dxi dxi dxi OXi 

dx' 3 dxi 

*3 — 5 -T-, = U 

0X1 dx 3 

dx[ dx 2 . dx 2 dx 2 

Wi = — - = — ■, = sin 0; m 2 =-—=-—, = cos 6; 

dx 2 dx[ dx 2 dx 2 

(144) 

dx 3 dx 2 
w 3 = = = 

OX 2 OX3 

dx[ dx 3 _ dx 2 dx 3 

«i = 3- = t-/ = 0; «2 =_- = __= 0; 

6^3 0*1 ox 3 ox 2 

dx 3 dx s 

ft* = T— = T1 r = I- 

da; 3 da;3 

Since ADP belongs to the orthorhombic bisphenoidal (Class 6), it will have 
the dielectric, piezoelectric and elastic tensors shown by equations (134), 
(135), (139). Applying equations (143) and (144) to these tensors it is 
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readily shown that the stresses for G = 45° are given by the equations ex- 
pressed in two index symbols 



rt = 



2 * 

+ (cfl + % 2 " 2C?6) Si + c? 3 s' t - hut 

(di + C12 — 2cu) r,' 

= _ ^i 



2 



, (cfi + ci 2 + 2c 66 ) / D / , , , (145) 

+ 2 T ^13>J8 -T «M03 

T"a = Ci3 5i + C13.S2 + CssSs 

T't = & S* + h u 5 2 ; £i = -huS't + 47r[/3n5i'] 

Tb = cu S' 6 - h u S[; £2 = hSx + 4irL8iifia] 

T't = (C " 7 gfz) S{ ; £3 = -*»[5j - 5j + 4»La»^. 

For a long thin longitudinally vibrating crystal all the stresses are zero 
except the stress T[ along the length of the crystal. Hence it is more ad- 
vantageous to use equations which express the strains in terms of the 
stresses since all the stresses can be set equal to zero except T x . All the 
strains are then dependent functions of the strain Si and this only has to 
be solved for. Furthermore, since plated crystals are usually used to 
determine the properties of crystals, and the field perpendicular to a plated 
surface is zero, the only field existing in a thin crystal will be £3 if the thick- 
ness is taken along the £3 or Z' axis. Hence the equations that express the 
strains in terms of the stresses and fields are more advantageous for calcu- 
lating the properties of longitudinally vibrating crystals. By 'orienting 
such crystals with respect to the crystallographic axis, all of the elastic 
constants except the shear elastic constants can be determined. All of 
the piezoelectric and dielectric constants can be determined from measure- 
ments on oriented longitudinally vibrating crystals. 

For such measurements it is necessary to determine the appropriate 
elastic, piezoelectric, and dielectric constants for a crystal oriented in any 
direction with respect to the crystallographic axes. We assume that the 
length lies along the x[ axis, the width along the x 2 axis and the thickness 
along the x'z axis. Starting with equations of the form 

Sij = SijktTkt + dij m Em 

t (146) 

5„= ^En + dntfTte 
4ir 
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and transforming to a rotated system of axes whose direction cosines are 
given by (142), the resulting equation becomes 

c' _ c * , dx '< dx 'i dx k dx * t'.j- j dXi dx 'i dXm *?' 
da?j da;,- da;* d^/ da;,- da;,- dXm 

s / e mn dx n dx m 1 ., dx n dXkdxt 1 

47T dx n OXm OX n OXk OX( 

All the stresses except T n can be set equal to zero and all the fields except 
£3 vanish. Furthermore, all the strains are dependently related to Sn . 
Hence for a thin longitudinal crystal the equation of motion becomes 

S' = s ?- /— — —— T' 4- d- dXl d* 1 dXm E' 
*' dxi dXi dx{ dx[ ,,m dXi dXi dxi m ' 

s r e mn dxz dx m _/ j_ , d^s da:* da;^ _,/ 
4?r dx„ d#8 da; n da;i da;i 

In terms of the two index symbols for the most general type of crystal, we 
have 

sun = su = sfid + (2sf 2 + saa)fiml + (2j? a + sfsKiWi 

+ 2(sf 4 + 5f 6 )^iWiWi + 2sf B £iWi + 2sf 6 /?Wi + sjf 2 w? 

+ (2*28 + su)mxni + 2s$ i miii 1 + 2(j? B + *fo)*»iA»i 

+ 25 2 S 6 Wi^i + sfs»i + 2S84»lWl + 2sf B «iA 

+ 2(5*6 + 5f B )«?Awi 
dm = dn = dnA^i + dvitwi + ^13^3^1 + dutgmini -f- dut3C1.n1 

+ duCjCinti + d2itn s Cl + J22W3W1 + d^-msnl + d^m^ini 

+ ditfitdtii + (h 6 m 3 Cmi + rf 3 iw 3 ^i + d 3 2n 3 mt + d^n^l 

+ dunstniUi + dstfhtith + d 36 n 3 Citni 

«33 = dm "f" 2tu*gmz "r* 2ei 3 ^3« 3 + €22^3 + 2e23W 3 » 3 + €33^3 

Hence by cutting 18 crystals with independent direction cosines 9 elastic 
constants and 6 relations between the remaining twelve constants can be 
determined. All of the piezoelectric constants and all of the dielectric 
constants can be determined from these measurements. These constants 
can be measured by measuring the resonant and antiresonant frequencies 
and the capacity at low frequencies. The resonant frequency f R is deter- 
mined by the formula 

11 ./I 



(149) 
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for any long thin crystal vibrating longitudinally. Hence when the density 
is known, sfl can be calculated from the resonant frequency and the length 
of the crystal. Using the values of s?i obtained for 15 independent orienta- 
tions enough data is available to solve for the constants of the first of 
equations (149). The capacities of the different crystal orientations meas- 
ured at low frequencies determine the dielectric constant e 33 and six orienta- 
tions are sufficient to determine the six independent dielectric constants 
tmn . The separation between resonance and antiresonance Af = f A — f R 
determines the piezoelectric constant d n according to the formula 



4. - \ a/M. i/? «ff . (i5i) 

The values of du measured for 18 independent orientations are sufficient 
to determine the eighteen independent piezoelectric constants. 

The remaining six elastic constants can be determined by measuring long 
thin crystals in a face shear mode of motion. Since this is a contour mode 
of motion, the equations are considerably more complicated than for a 
longitudinal mode and involve elastic constants that are not constant field 
or constant displacement constants. It can be shown that the fundamental 
frequency of a crystal with its length along X\ , width (frequency determining 
direction) along xo and thickness (direction of applied field) along X3 , will be 

1 / C,E I C,E , - / 1 C,E C,B\2 1 * C.E* 

f = — A/ ° 22 ~L ° 66 ^ ggg ~ 66 I Z 26 (152) 

2( w » 2p 

where the contour elastic constants are given in terms of the fundamental 
elastic constants by 



(153) 
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(154) 



Since all of the constants except $f 2 and $f 6 can be determined by measure- 
ments on longitudinal crystals and the value of (2$f 2 + $6e) has been de- 

6 This is proved in a recent paper "Properties of Dipotassiura Tartrate (DKT) Crys- 
tals," Phys. Rev., Nov., 1946. 
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termined, the measurement of the lowest mode of the face shear crystal 
gives one more relation and hence the values of sf 2 and s% 6 can be uniquely 
determined. 

Similar measurements with crystals cut normal to xi and width along #3 
and with crystals cut normal to #2 and width along X\ determine the constants 
Su , Sm and sft , Sn respectively. The equivalent constants are obtained 
by adding 1 to each subscript 1, 2, 3 or 4, 5, 6 for the first crystal with the 
understanding that 3+1=1 and 6+1 = 4. For the second crystal 2 
is added to each subscript. 

Finally the remaining three constants can be determined by measuring 
the face shear mode of three crystals that have their lengths along one of 
the crystallographic axes and their width (frequency determining axis) 
45° from the other two axes. 

Any symmetry existing in the crystal will cut down on the number of 
constants and hence on the number of orientations to determine the funda- 
mental constants. 

6. Temperature Effects in Crystals 

In section 2 a general expression was developed for the effects of tempera - 
ture and entropy on the constants of a crystal. Two methods were given, 
one which considers the stresses, field, and temperature differentials as the 
independent variables, and the second which considers the strains, displace- 
ments and entropy as the independent variables. In tensor form the 10 
equations for the first method take the form 

T k C = c%USij - h" nk ib n - X?/ dQ 
E m =- h' m 1 jS i j + 47T/&; 8 n - & D dQ (155) 

dQ = -Q\% D Sa - eq* D 8 n + "^ 

The piezoelectric relations have already been discussed for adiabatic condi- 
tions assuming that no increments of heat dQ have been added to the 
crystal. 

If now an increment of heat dQ is suddenly added to any element of the 
crystal, the first equation shows that a sudden expansive stress is generated 

S D 

proportional to the constant \l'( which has to be balanced by a negative 
stress (a compression) in order that no strain or electric displacement shall 
be generated. This effect can be called the stress caloric effect. The 
second equation of (155) shows that if an increment of heat dQ is added to 
the crystal an inverse field E m has to be added if the strain and surface 
charge are to remain unchanged. This effect may be called the field caloric 
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effect. Finally the third equation of (155) shows that there is a reciprocal 
effect in which a stress or a displacement generates a change in temperature 
even in the absence of added heat dQ. These effects can be called the strain 
temperature and charge temperature effects. 

The second form of the piezoelectric equations given by (58) are more 
familiar. In tensor form these can be written 

Sii = Si/UTtf + d%ijE m + af,- dQ 
8„ = G*tTU + ^ E m + pi dQ (156) 

d Q = Qd a= Q a E k( T k l + OplEm + pC* p dQ 

The a-fy are Ihe temperature expansion coefficients measured at constant 
field. In general these are a tensor of the seond rank having six com- 
ponents. The constants p T n are the pyroelectric constants measured at 
displacements which relate the increase in polarization or surface charge 
due to an increase in temperature. They are equal to the so-called "true" 
pyroelectric constants which are the polarizations at constant volume caused 
by an increase in tempeiature plus the "false" pyioelectric effect of the 
first kind which represents the polarization caused by a uniform temperature 
expansion of the crystal as its temperature increases by dQ. As mentioned 
previously it is more logical to call the two effects the pyroelectric effects 
at constant stress and constant strain. By eliminating the stresses from 
the first of equations (156) and substituting in the second equation it is 
readily shown that 

pl = pi - a'Aa ( 157 ) 

Hence the difference between the pyroelectric effect at constant stress and 
the pyroelectric effect at constant strain is the so-called "false" pyroelectric 
effect of the first kind «»,-«»»,• . 

The first term on the right side of the last equation is called the heat of 
deformation, for it represents the heat generated by the application of the 
stresses Tul ■ The second term is called the electrocaloric effect and it 
represents the heat generated by the application of a field. The last term 
is p times the specific heat at constant pressure and constant field. 

The temperature expansion coefficients a,-, form a tensor of the second 
rank and hence have the same components for the various crystal classes 
as do the dielectric constants shown by equation (134). 

The pyroelectric tensor p T n and p s „ are tensors of the first rank and in 
general will have three components pi, Pz, and pi . In a similar manner 
to that used for second, third and fourth rank tensors it can be shown that 
the various crystal classes have the following components for the first rank 
tensor p„ . 
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Class 1 : components pi , p2 , pa • 

Class 3 : Y axis of binary symmetry, components 0, p 2 , (158) 

Class 4: components pi , 0, p 3 . 

Classes 7, 10, 14, 16, 20, 23, and 26: components 0, 0, p 3 ; and Classes 
2, 5, 6, 8, 9, 11, 12, 13, 15, 17, 18, 19, 21, 22, 24, 25, 27, 28, 29, 30, 31, and 
32 : components 0, 0, 0, i.e., p = 0. 

For a hydrostatic pressure, the stress tensor has the components 

T n = T 22 = r 33 = -p = pressure; T u = r 13 = 7^ = (159) 
Hence the displacement equations of (156) can be written in the form 

8 n = e i^M m -d e n p + p T n dQ (160) 

4tt 

where 

Onp = dniiTn + ^22^22 + OutsT^ 

that is the contracted tensor d n kkTkk • This is a tensor of the first rank 
which has the same components as the pyroelectric tensor p n for the various 
crystal classes. 

7. Second Order Effects est Piezoelectric Crystals 
We have so far considered only the conditions for which the stresses and 
fields are linear functions of the strains and electric displacements. A 
number of second order effects exist when we consider that the relations are 
not linear. Such relations are of some interest in ferroelectric crystals such 
as Rochelle salt. A ferroelectric crystal is one in which a spontaneous 
polarization exists over certain temperature ranges due to a cooperative 
effect in the crystal which lines up all of the elementary dipoles in a given 
"domain" all in one direction. Since a spontaneous polarization occurs in 
such crystals it is more advantageous to develop the equations in terms of 
the electric displacement rather than the external field. Also heat effects 
are not prominent in second order effects so that we develop the strains and 
potentials in terms of the stresses and electric displacements D. By means 
of McLaurin's theorem the first and second order terms are in tensoi form 

Sii = dTd U + d8 n 8n + 2! \_dT k ldT qr TktTqr 

5„a5o M °J 



d 2 S-- d 2 S- 

+ 2 ^F7t r *^» + ^717 5 » 5 ° I + • • * ni g her terms 

8E m dE m 1|~ d 2 E m 

*" " dTkC Tkt + d8 n 8n + 2! IdTddT* TktTqr 

+ 2 S£k ™* + S^o 8J °] + • • • higher terms 

whereas before 5 = D/4ir. 



(161) 



• 
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In this equation the linear partial differentials have already been discussed 
and are given by the equations 

dTkl - s ijkt ; ^ - dT __ - g iin , d8n - *irp mn UO/J 

where sljU ar e the elastic compliances of the crystal at constant displace- 
ment, gn n the piezoelectric constants relating strain to electric displacement 
/4t, and (3m n the dielectric "impermeability" tensor measured at constant 
stress. We designate the partial derivatives 

ftgg _ „i> , ggg = ^E m ^ 

dTkidT qr iyiikt " r ' dT k (dd m dT k 4dT qT Miiktn 

(163) 
d Sij _ d E n __ „o d E m _qD 

dd n d8o BTijdSo ' B8 n d8o 

The tensors N, M, Q, and are respectively tensors of rank 6, 5, 4 and 3 
whose interpretation is discussed below. Introducing these definitions 
equations (161) can be written in the form 

Sij= T k ({ S D ijk t + \N D iik t qT T qr + M ° ik (n8n] + Bnlgijn + Wlnoh] 

(164) 

Em = T k t\g mk l+\M D HklnT qT + Qk^nSn) + S„[4r^ n + \OLoh] 

Written in this form the interpretation of the second order terms is obvious. 
Nljktqr represents the nonlinear changes in the elastic compliances $?,*/ 
caused by the application of stress T qr . Since the product of Niikl qT T qr 
represents a contracted fourth rank tensor, there is a correction term for 
each elastic compliance. The tensor M^ jk e n can represent either the non- 
linear correction to the elastic compliances due to an applied electric dis- 
placement D„ or it can represent the correction to the piezoelectric constant 
gij„ due to the stresses Trf . By virtue of the second equation of (162), 
the second equivalence of (163) results. The fourth rank tensor \Qij no 
represents the electrostrictive effect in a crystal for it determines the strains 
existing in a crystal which are proportional to the square of the electric 
displacement. Twice the value of the electrostrictive tensor \Q°j no , which 
appears in the second equation of (164) can be interpreted as the change 
in the inverse dielectric constant or "impermeability" constant. Since a 
change in dielectric constant with applied stress causes a double refraction 
of fight through the crystal, this term is the source of the piezo-optical effect 
in crystals. The third rank tensor 0m„ o represents the change in the "im- 
permeability" constant due to an electric field and hence is the source of 
the electro-optical effect in crystals. 

These equations can also be used to discuss the changes that occur in 
ferroelectric type crystals such as Rochelle Salt when a spontaneous polariza- 
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tion occurs in the crystal. When spontaneous polarization occurs, the 
dipoles of the crystal are lined up in one direction in a given domain. For 
Rochelle salt this direction is the ±.Y axis of the crystal. Now the electric 
displacement D x is equal to 

where P Xo is the electronic and atomic polarization, and P XD the dipole 
polarization The electronic and atomic polarization is determined by the 
field and hence can be combined with the field through the dielectric constant 
to , which is the temperature independent part of the dielectric constant. 
When the crystal becomes spontaneously polarized, a field E x will result, but 
this soon is neutralized by the flow of electrons through the surface and 
volume conductance of the crystal and in a short time E x = 0. Hence for 
any permanent changes occurring in the crystal we can set 

8 X = — - =Px D = dipole polarization (166) 

47T 

which we will write hereafter as P\ . 

In the absence of external stresses the direct effects of spontaneous polari- 
zation are a spontaneous set of strains introduced by the product of the 
spontaneous polarization by the piezoelectric constant, and another set 
produced by the square of the polarization times the appropriate electro- 
strictive components. For example, Rochelle salt has a spontaneous 
polarization P\ along the X\ axis between the temperatures — 18°C to 
+ 24°C. The curve for the spontaneous polarization as a function of 
temperature is shown by Fig. 6. The only piezoelectric constant causing 
a spontaneous strain will be gu/t = gm ■ Hence the spontaneous polariza- 
tion causes a spontaneous shearing strain 

S, = gu P x = 120 X 10~ 8 X 760 = 9.1 X 10" 4 (W 

if we introduce the experimentally determined values. Since S* is the 
cosine of 90° plus the angle of distortion, this would indicate that the right 
angled axes of a rhombic system would be distorted 3.1 minutes of arc. 
This is the value that should hold for any domain. For a crystal with 
several domains, the resulting distortion will be partly annulled by the 
different signs of the polarization and should be smaller. Mueller measured 
an angle of 3'45" at 0°C for one crystal. This question has also been 

6 This has been measured by measuring the remanent polarization, when all the domains 
are lined up. See "The Dielectric Anomalies of Rochelle Salt," H. Mueller, Annals of 
tlu: N. Y. Acad. Science, Vol. XL, Art. 5, page 338, Dec. 31, 1940. 

7 "Properties of Rochelle Salt," H. Mueller, Phys. Rev., Vol. 57, No. 9, May 1, 1940. 
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investigated by the writer and Miss E. J. Armstrong by measuring the 
temperature expansion coefficients of the Y and Z axes and comparing their 
average with the expansion coefficient at 45° from these two axes. The 
difference between these two expansion coefficients measures the change 
in angle between the Y and Z axes caused by the spontaneous shearing 
strains. The results are shown by Fig. 7. Above and below the ferro- 
electric region, the expansion of the 45° crystal coincides with the average 
expansion of the Y and Z axes measured from 25°C as a reference tempera- 
ture. Between the Curie temperatures a difference occurs indicating that 
the Y and Z crystallographic axes are no longer at right angles. The dif- 
ference in expansion per unit length at 0°C (the maximum point) corresponds 
to 1.6 X 10 -4 cm per cm. This represents an axis d istortion of 1.1 minutes 
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Fig. 6. — Spontaneous polarization in Rochelle Salt along the X axis. 

of arc. Correspondingly smaller values are found at other temperatures 
in agreement with the smaller spontaneous polarization at other tempera- 
tures. It was also found that practically the same curve resulted for either 
45° axis, indicating that the mechanical bias put on by the optometer used 
for measuring expansions introduced a bias determining the direction of the 
largest number of domains. 

The second order terms caused by the square of the spontaneous polariza- 
tion is given by the expression 

Sa = QimPl (168) 

Since Q is a fourth rank tensor the possible terms for an orthorhombic 
bisphenoidal crystal (the class for Rochelle salt) are 

Sn = QnnPl J Sn - (guA ; S„ = QnnP\ (169) 
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In an effort to measure these effects, careful measurements have been made 
of the temperature expansions of the three axes X, Y and Z. The results 
are shown by Table II. On account of the small change in dimension from 
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Fig. 7. — Temperature expansion curve along an axis 45° between Y and Z as a 
function of temperature. 



the standard curve it is difficult to pick out the spontaneous components 
by plotting a curve. By expressing the expansion in the form of the 
equation 



^ = ai(T-2S) + a^T-lSf + a 3 (T-25)' 



(170) 
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Table II 
Measured Temperature Expansions for the Three Crystalographic Axes 



Temperature Expansion 


Temperature 
in °C. 


Expansion 

x io-« 

Y Axis 


Temperature 
in°C. 


Expansion 


in °C. 


X H>-' 
AT Axis 


x io-« 

Z Axis 


39.6 
38.7 
35.2 


10.2 
9.46 
6.96 


+35.0 
30.3 
25.25 


4.45 

2.5 

0.2 


+34.9 
29.9 
25.05 


+4.9 
2.5 

+ .05 


30.2 
27.2 
26.2 


3.63 
1.41 
0.71 


23.9 
22.9 
19.35 


-0.42 
-0.88 
-2.4 


24.0 

19.95 

14.95 


-.5 
-2.62 
-5.11 


25.15 

24.0 

23.0 


0.06 
-0.71 
-1.39 


14.9 

10.0 

5.4 


-4.25 
-6.25 
-8.18 


+9.75 

+4.9 




-7.55 

-9.9 

-12.31 


21.8 
16.0 
15.2 


-2.37 

-6.5 

-7.05 


+0.3 

-9.7 

-16.3 


-10.15 
-13.98 
-16.41 


-6.35 

-10.5 
-15.0 


-15.3 

-17.29 

-19.42 


4.9 
+0.3 
-4.7 


-14.12 
-17.28 
-20.3 


-20.85 

-25.1 

-30.3 


-17.94 
-19.22 
—20.8 


-18.0 
-23.2 
-25.1 


-20.86 
-23.08 
-23.96 


-10.7 
-15.3 
-20.7 


-24.0 
-26.6 
-30.2 


-35.0 
-39.7 

-53.2 


-22.23 
-23.54 
-27.60 


-31.1 
-35.0 
-40.0 


-26.59 
-28.28 
-30.4 


-25.7 
-30.1 
-34.7 


-32.7 
-35.2 
-37.85 










-40.7 
-45.0 
-50.5 


-41.25 

-44.0 

-47.0 











and evaluating the constants by employing temperatures outside of the 
ferroelectric range, a normal curve was established. For the X, Y, and Z 
axes these relations are 

^ = 69.6 X l(T 6 (r-25) + 7.4 X l(T 8 (T-25) 2 - 3.13 X l(T 10 (T-25) 3 

■Lt 

(X direction) 



— = 43.7 X l(T 6 (T-25) + 8.16 X l(T 8 (T-25) 2 .- 3.60 X 10 10 (T-25) 3 

(F direction) 



(171) 



$£. = 49.4 X 10" 6 (r-25) + 1.555 X 10" 8 (r-25) 2 - 2.34 X 10- 10 (r-25) 3 

(Z direction) 
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The difference between the normal curves and the measured values in the 
Curie region is shown plotted by the points of Fig. 8. The solid and dashed 
curves represent curves proportional to the square of the spontaneous 
polarization and with multiplying constants adjusted to give the best fits 
for the measured points. These give values of Qitu , Qaau , Qsm equal to 

Qim = -86.5 X 1(T 12 ; Q? 2U = +17.3 X 10"*; 

Qnn= -24.2 X1(T 12 (172) 

Another effect noted for Rochelle salt is that some of the elastic constants 
suddenly change by small amounts at the Curie temperatures. This is a 
consequence of the tensor M?,^„ , for if a spontaneous polarization P 
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Fig. 8.— Spontaneous electrostrictive strain in Rochelle Salt along the 
three crystallographic axes. 

occurs, a sudden change occurs in some of the elastic constants as can be seen 
fiom the first of equations (164). The second equation of (164) shows 
that this same tensor causes a nonlinear response in the piezoelectric con- 
stant. Since a change in the elastic constant is much more easily deter- 
mined than a nonlinear change in the piezoelectric constant, the first effect 
is the only one definitely determined experimentally. Since all three crys- 
tallographic axes are binary axes in Rochelle salt, it is easily shown that 
the only terms that can exist for a fifth rank tensor are terms of the types 

Mum ; Mf 2 o 2 3 ; Mf 23 33 (173) 

with permutations and combinations of the indices. Hence when a spon- 
taneous polarization Pi occurs, the elastic constants become 

s D iik t - IffrftPi (174) 
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Comparing these with the relation of (90) we see that the spontaneous 
polarization has added the elastic constants 

d (Af 11231 + Af 11321 + Af 23111 + Af32 in)P i 
5u= 2 

v _ (Af ?2231 + Af 22321 + Af 23221 + Af 3 222l)-Pl 
$24 — n " 



d (M£j33i + M%m\ + A/33321 + MtzmjPi 

$34 = ~ 



(175) 



(Ml21i.l + Af 13211 + Af 31121 + Af 3 1211 
D + Aff 2 131 + Af 12311 + Af 21131 + M 2131l)-Pl 

$66 — 2 

between the two Curie points. Hence while the spontaneous polarization 
Pi exists, the resulting elastic constants are 



$11 , 


$12 , 


$13 , 


$14 , 


, 





Sn , 


$22, 


$23 , 


$2 J , 


, 





Su, 


$23, 


$33 , 


$34 , 


, 





$14 , 


$24 , 


$34 , 


$44 , 


, 





, 
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$55 , 


$56 
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, 


, 


, 


$56 , 


$C6 



(176) 



Comparing this to equation (139) which shows the possible elastic constants 
for the various crystal classes, we see that between the two Curie points, 
the crystal is equivalent to a monoclinic sphenoidal crystal (Class 3) with 
the X axis the binary axis. Outside the Curie region the crystal becomes 
orthorhombic bisphenoidal. This interpretatiorTagrees with the tempera- 
ture expansion curves of Fig. 7. 

The sudden appearance of the polarization Pi will affect the frequency 
of a 45° X-cut crystal, for with a crystal cut normal to the X axis and with 
the length of the crystal at an angle with the Y axis, the value of the 
elastic compliance $22 along the length is 
4' = 4 cos 4 6 + 2su cos 3 sin 6 + (2sn + sfi) sin 2 cos 2 

+ 2$f 4 sin 3 cos + $£, sin 4 

Hence for a crystal with its length 45° between the Y and Z axes, elastic 
compliance becomes 



D 1 r,/ D 1 D I »\ 1 Jt \ D 

ID _ $22 + 2 ($24 + $2 3 + $34 ) + $44 + $33 
$22 — ~. 



(178) 
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For a 45° X-cut crystal we would expect a sudden change in the value of 
sn as the crystal becomes spontaneously polarized between the two Curie 
points due to the addition of the 524 and su elastic compliances. Such a 
change has been observed for Rochelle salt 8 as shown by Fig. 9 which shows 
the frequency constant of a nonplated crystal for which the elastic com- 
pliances sfj should hold. 
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Fig. 9.— Frequency constant and Q of an unplated 45°X cut Rochelle Salt crystal 
plotted as a function of temperature. 

Hence the sudden change in the elastic constant is a result of the two 
second order terms J& + 5? 4 , which are caused by the spontaneous polariza- 
tion. The value of the sum of these two terms at the mean temperature 
of the Curie range, 3°C is 

5 ° + ,£ = 4.1 X 1(T M cmVdyne (179) 

Crystals cut normal to the Y and Z axes should not show a spontaneous 
change in their frequency characteristic since the spontaneous terms s u , 
s u , 5 34 and 5 66 do not affect the value of Young's modulii in planes normal 
to Y and Z. Experiments on a 45° F-cut Rochelle salt crystal do not show 
a spontaneous change in frequency at the Curie temperature, although there 
is a large change in the temperature coefficient of the elastic compliance 
between the two Curie points. This is the result of third order term and is 

8 "The Location of Hysteresis Phenomena in Rochelle Salt Crystals," W. P. Mason, 
Phys. Rev., Vol. 50, p. 744-750, October 15, 1940. 
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not considered here. The spontaneous s 66 constant affects the shear con- 
stant s 66 for crystals rotated about the X axis and could be detected experi- 
mentally. No experimental values have been obtained. 

The effects of spontaneous polarization in the second equation of (164) 
are of two sorts. For an unplated crystal, a spontaneous voltage is gen- 
erated on the surface, which, however, quickly leaks off due to the surface 
and volume leakage of the crystal. The other effects are that the spon- 
taneous polarization introduces new piezoelectric constants through the 
tensor Qk/m,, , changes the dielectric constants through the tensor 0™„ and 
introduces a stress effect on the piezoelectric constants through the tensor 
M k f m gr . Since piezoelectric constants are not as accurately measured as 
elastic constants, the first effect has not been observed. The additional 
piezoelectric constants introduced by the tensor Qkpnn are shown by equa- 
tion (180) 

gu gu gta gu 

g 2b g 26 (180) 

£35 g30 

Since the only constants for the Rochelle salt class, the orthorhombic 
bisphenoidal, are gu , g 2 6 , gs6 , this shows that between the two Curie points 
the crystal becomes monoclinic sphenoidal, with the X axis being the 
binary axis. The added constants are, however, so small that the accuracy 
of measurement is not sufficient to evaluate them. From the expansion 
measurements of equation (172) and the spontaneous polarization values, 
three of them should have maximum values of 

g n = -6.6 X HT 8 ; g u - +1.3 X 10" 8 ; gu = -1.8 X 10 -8 (181) 

These amount to only 6 per cent of the constant gu , and hence they are 
not easily evaluated from piezoelectric measurements. 

The effect of the tensor Om„„ is to introduce a spontaneous dielectric 
constant €23 between the Curie temperatures so that the dielectric tensor 
becomes 

eu, , 

0, C22, 623 (182) 

, €23 , 633 

As discussed at length by Mueller 9 this introduces a spontaneous bire- 
fringence for light passing through the crystal along the X, Y and Z axes 
which adds to the birefringence already present. 

9 "Properties of Rochelle Salt I and IV," Phys. Rev. 47, 175 (1935) ; 58, 805 November 1, 
1940. 



